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ABSTRACT 


The  ordinary  diffusion  coefficient  and  the  Soret  (thermal 
diffusion)  coefficient  were  measured  for  a  solution  of  3.5  weight  percent 
sodium  chloride  in  water;  this  system  was  selected  as  a  reasonable  binary 
simulant  to  the  multicomponent  sea  water  system.  Ordinary  diffusion  co¬ 
efficients  were  measured,  at  25°C,  for  pressures  over  the  range  1  to  1000 
bars;  the  Soret  coefficient  was  determined  at  27°C  and  approximately 
1000  bars  pressure,  and  compared  to  measured  values  at  1  bar  pressure. 

The  experimental  high  pressure  cells  and  techniques  developed  for  measuring 
these  properties  for  binary  salt  solutions  can  be  modified  and  adapted  for 
possible  analogous  measurements  with  multicomponent  solutions. 

The  influence  of  pressure  on  the  ordinary  diffusion  coefficient 
in  3.5  weight  percent  NaCl-I^O  at  25°C  was  found  to  be  negligible  over 
the  pressure  range  1  to  lOOCTbars.  In  view  of  the  widely  differing 
diffusional  behavior  of  the  various  binary  salt  constituents,  however, 
a  significant  pressure  effect  could  exist  for  a  multicomponent  mixture  such  as  sea 
water.  Based  upon  the  data  from  low  pressure  experiments  and  a  single 
high  pressure  experiment  in  the  Soret  diffusion  apparatus,  the  Soret  co¬ 
efficient  at  room  temperature  likewise  appears  to  be  virtually  independent 
of  pressure  up  to  1000  bars.  The  extrapolated  reversal  of  direction  for 
thermal  diffusion  in  the  NaCl-H-O  system  at  temperatures  lower  than  about 
7°C  should,  however,  be  investigated  at  high  pressures  for  possible  signi¬ 
ficance  regarding  the  sea  water  system  itself. 

An  extensive  review  and  analysis  was  made  of  the  available  lit¬ 
erature  pertaining  to  the  viscosity  and  thermal  conductivity  of  multicomponent 
salt  solutions,  including  sea  water.  This  analysis  has  shown  that  these 
properties  are  not  too  different  from  those  of  pure  water,  and  can  be  pre¬ 
dicted  satisfactorily  for  engineering  purposes  by  simple  additivity 
relationships  based  upon  the  most  important  salt  ions  present.  Neither 
the  viscosity  nor  the  thermal  conductivity  appears  to  exhibit  any  unusual 
effect  (relative  to  water)  with  changes  in  temperature,  pressure,  or  con¬ 
centration,  and  no  direct  relationships  apparently  exist  between  either 
viscosity  or  thermal  conductivity  and  the  mass  diffusional  properties  of 
the  individual  ions. 
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SECTION  1 


INTRODUCTION 


1.1  PRELIMINARY  CONSIDERATIONS 

There  are  a  number  of  situations  of  interest  to  physical  and 

chemical  oceanography  where  the  molecular  diffusion  of  ions  in  an  aqueous 

solution  can  be  important.  These  situations  include,  for  example,  the 

molecular  transport  phenomena  that  occur  in  the  vicinity  of  the  interface 

between  sea  water  and  sea  ice,  between  undersea  instruments  or  vessels 
* 

and  sea  water  ,  and  between  the  ocean  bottom  and  sediment  regions  and  the 
surrounding  sea  water.  Important  technological  problems  such  as  the  rate 
of  corrosion  of  metals  by  sea  water  and  the  performance  of  undersea 
batteries  are  also  significantly  affected  by  molecular  ionic  diffusion 
processes. 

The  two  most  common  causes  for  ion  diffusion  in  these  situations 
are  concentration  gradients  and  temperature  gradients.  The  ion  flux  is 
directly  related  to  these  two  driving  forces  by  molecular  diffusion 
coefficients  known  as  the  ordinary  diffusion  coefficient  and  the 


*  In  the  laminar  boundary  layer,  or  in  the  viscous  sublayer  of  a  turbulent 
boundary  layer  where  molecular  processes  are  of  significance. 

**  Often  called  isothermal 
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Soret  coefficient.  These  coefficients,  in  general,  depend  upon  con¬ 
centration,  temperature,  and  pressure.  Their  definition  for  binary  liquid 
systems  (for  example,  salt  +  water),  whether  concentrated  or  dilute,  has 
been  well  established  (see  Section  1.3).  The  theoretical  prediction  of 
these  coefficients  is  only  possible  for  some  dilute  binary  systems  using 
previously  determined  property  parameters  which  are  usually  measured  at 
temperatures  near  25°C  and  at  a  pressure  of  one  atmosphere.  Generally ,  pre¬ 
dictions  for  concentrated  binary  solutions  and  for  multicomponent 
solutions  are  >not  now  possible.  Therefore,  it  is  necessary  to  measure 
these  coefficients  under  well  defined  conditions. 

The  molecular  diffusion  properties  of  sea  water,  a  somewhat 
concentrated  multicomponent  salt  solution,  are  not  known.  The  binary 
aqueous  solution  diffusion  properties  of  its  major  salt  constituent, 
sodium  chloride,  are  known  only  for  temperatures  near  25°C  and  for  a 
pressure  of  1  atmosphere.  Evan  less  is  known  about  these  properties  for 
the  other  salt  constituents  of  sea  water. 

Since  temperatures  in  the  oceans  can  vary  from  -2°C  to  30°C,  and 
since  pressure  can  vary  with  depth  from  1  atmosphere  to  approximately  1000 
atmospheres,  the  diffusion  properties  should  be  known  within  these  limits 
of  temperature  and  pressure. 

It  is  known  that  the  diffusion  properties  of  multicomponent  salt 
solutions  can  be  much  different  than  those  of  binary  salt  solutions.  Both 
the  direction  of  diffusion  and  the  diffusion  rates  can  be  altered  depending 
upon  the  relative  concentrations  of  the  various  types  of  lone  in  the  multi- 
component  system.  Clearly,  a  detailed  understanding  of  molecular  diffusion 
phenomena  in  such  complex  systems  can  only  be  obtained  after  a  reasonably 
good  description  of  diffusion  in  the  various  component  binary  systems  has 
been  developed. 

*  Named  after  an  early  investigator  of  thermal  diffusion  phenomena  in 
aqueous  salt  solutions. 
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The  following  sections  first  briefly  outline  the  objectives 
of  the  research  investigation  reported  herein  and  give  the  equations 
which  exactly  define  the  diffusion  properties  of  interest  in  this 
study. 

1.2  OBJECTIVES  OF  THIS  INVESTIGATION 

In  view  of  the  discussion  of  the  preceding  section,  it  is 
appropriate  to  consider  the  diffusion  properties  of  binary  systems  first, 
and  in  particular  those  of  the  sodium  chloride— water  system,  before 
attempting  to  study  the  more  complex  case  of  sea  water  itself.  Because 
sodium  and  chloride  ions  are  the  most  predominate  ions  in  sea  water,  a 
3.5  percent  by  weight  solution  of  NaCl-H^O  is  a  convenient  first 
approximation  to  typical  sea  water  (7). 

On  this  basis,  the  following  objectives  were  selected  for 
this  investigation: 

(i)  to  measure  the  binary  Soret  coefficient  and  the  ordinary  diffusion 
coefficient  of  the  synthetic  sea  water,  3.5  percent  by  weight 
sodium  chloride  in  water,  and  to  correlate  the  results  with 
known  theoretical  considerations; 

(ii)  to  investigate  the  ordinary  and  thermal  diffusion  properties 

of  binary  aqueous  salt  solutions  other  than  the  NaCl-^O  system; 

(iii)  to  develop  an  experimental  apparatus  for  measuring  the  high 
pressure  thermal  diffusion  properties  of  multicomponent  salt 
solutions ; 

(iv)  to  determine  how  the  experimental  apparatus  used  to  measure 
the  binary  ordinary  diffusion  coefficient  can  be  modified  in 
order  to  measure  the  ordinary  diffusion  properties  of  multi* 
component  aqueous  electrolyte  solutions; 

(v)  to  relate  analytical  predictions  to  experimental  results; 


*  Numbers  in  parentheses  which  occur  in  the  text  refer  to  the  references 
listed  at  the  end  of  this  report. 
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(vi)  to  review  and  analyze  the  available  literature  pertaining  to 
the  viscosity  and  thermal  conductivity  of  multicomponent 
aqueous  salt  solutions,  including  sea  water,  and  to  relate, 
where  applicable,  this  property  information  to  the  mass 
diffusion  property  data. 

The  scope  of  the  investigation  defined  by  these  objectives 
does  not  include  measurements  of  the  diffusion  properties  of  multi- 
component  systems.  An  important  aspect  of  the  study  will  be  to  determine 
how  pressure  affects  the  ordinary  diffusion  coefficient  and  Soret 
coefficient  of  the  sodium  chloride— water  system. 


1.3  PHENOMENOLOGICAL  DESCRIPTION 

The  diffusion  coefficients  are  defined  by  general  flux-driving 
force  equations  of  the  following  form  (for  binary  systems — gaseous,  two 
organic  liquids,  or  one  salt  in  an  electrically  neutral  solvent,  VP  -  0, 
and  no  differing  body  forces  acting  except  for  the  interionic  electrical 
forces)  (49;  Eqs.  3.26  and  3.27)  (50): 


j/  *  *CDi2  {7x2  +  CT12X1X2VT} 

or, 

V  *  -pD12  {vy2  +  C7i2y1y2vT} 


(1) 


(2) 


where  D^2  denotes  the  binary  ordinary  diffusion  coefficient  and  a ^  denotes 
the  binary  Soret  coefficient  (symbols  are  defined  in  the  nomenclature 
section).  The  fluxes  are  related  to  species  velocities  by  the  equations: 


(3) 
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Equations  1  and  2  define  the  same  diffusion  coefficients  and  can  be  shown 
to  be  identical. 

The  so-called  Fick  diffusion  coefficient  is  defined  by  Pick's 
first  law  (isothermal  conditions) 


'  °Fick 


(7) 


where , 


C2V2 


(8) 


The  inadequacy  of  an  ordinary  diffusion  coefficient  defined  in  this  manner 
and  its  nonequivalence  with  have  been  thoroughly  discussed  in  the 
literature  (51). 

A  primary  objective  of  this  investigation  is  to  learn  more  about 
the  pressure  and  temperature  dependence  of  and  <j  ry  for  the  3.5  weight 
percent  sodium  chloride— water  system. 

The  basic  defining  relationships  take  on  somewhat  different  forms 
for  multicomponent  ionic  systems.  These  will  be  discussed  in  the  sections 
that  follow. 
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SECTION  2 


ORDINARY  DIFFUSION  COEFFICIENT  STUDY 


2 . 1  ANALYTICAL  CONSIDERATIONS 

The  most  elementary  description  of  ordinary  diffusion  in 
liquids,  is  given  by  the  Stokes-Einstein  relation  (50,  pp.  513-515), 


RT 


12  6tt  m.,R0N 
2  o 


(9) 


for  the  diffusion  of  spherical  particles  of  such  a  size  that  the  solvent 
(species  1)  appears  to  the  diffusing  species  as  a  continuum.  This 
relation  predicts  a  temperature  dependence  of 


Dj2  ~(T/m.j) 


(10) 


which  is  borne  out  by  experimental  results  for  certain  systems,  such  as 
the  dilute  NaCl-I^O  system. 

For  very  dilute  (x2~*o),  binary  electrolyte  solutions,  the 
Nernst  expression, which  neglects  interionic  forces,  has  some  utility. 

This  simple  relation  relates  D^  (D^  f°r  x2“*  °)  t0  the  limiting  equivalent 
conductivities  by  the  equation  (49), 
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A  more  elaborate  theory  has  been  developed  by  Onsager  and 
Fuoss  (52}  pp. 105-122,  243-255)  and  is  more  generally  useful  for 
moderately  dilute  systems.  This  theory  is  inadequate  for  concentrated 
systems,  as  will  be  seen  below. 

Gordon  (53)  has  proposed  a  semi-empirical  expression  in  order 
to  extrapolate  In  the  very  dilute  range  to  high  concentrations. 

Very  little  use  has  been  made  of  Gordon's  relation  even  though  quite 
satisfactory  agreement  between  values  calculated  from  it  for  high  con¬ 
centrations  and  experimentally  measured  values  has  been  obtained.  Van 
Rysselberghe  (54)  has  proposed  a  similar  expression. 

The  temperature  dependence  of  for  the  NaCl-^O  system  at 
infinite  dilution  is  shown  in  Figure  2.1.  The  D°2  values  were  calculated 
from  Eq.  11  using  the  limiting  equivalent  conductivities,  X°,  in  (56,  p.  465). 
The  viscosities  used  for  the  abscissa  are  values  reported  for  pure  water 
(56,  p.457).  The  temperature  dependence  predicted  by  Eq.  10  is  followed 
nearly  exactly.  It  is  apparent  that  the  temperature  dependence  of  Dj^  Is 
significant.  For  the  case  given  in  Figure  2.1,  there  is  a  three-fold 
increase  in  in  going  from  5°C  to  55°C. 

All  published,  reliable,  ordinary  diffusion  coefficient  data 
for  the  NaCl-I^O,  MgC^-l^O,  and  CaC^-HL^O  systems  as  a  function  of  composition 
and  temperature  at  1  atm  is  presented  in  Figures  2.2,  2.3  and  2.4.*  A  smooth 

*  The  sources  for  the  data  are  ns  follows: 

D  Q_Narl :  Harned  &  Hildreth  (57),  Vitagliano  &  Lyons  (58),  Dunlop  &  Gosting 
H2  ’  (59),  Stokes  (60,  O'Donnell  &  Gosting  (61),  Gordon  (53),  Clack  (62) 


D„  n  «i  :  Harned  &  Polestra  (63),  Oholm  (64) 

H20-MgCl2 

DH  0-CaCl  :  Harned  &  tevy  (65);  Hall,  Wishaw,  &  Stokes  (66);  Lyons  &  Riley  (67) 
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FIGURE  2.1.  TEMPERATURE  DEPENDENCE  OF  THE  ORDINARY 
DIFFUSION  COEFFICIENT  FOR  THE  NaCl-H20 
SYSTEM  AT  INFINITE  DILUTION  (CALCULATED 
FROM  EXPERIMENTALLY  DETERMINED  LIMITING 
EQUIVALENT  CONDUCTIVITIES) (1  ATM) 
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NERNST 


ORDINARY  DIFFUSION  COEFFICIENT  FOR  THE  NaCl-I^O 


NERNST  LIMITING 
VALUE  (25°C) 


FIGURE  2.3.  ORDINARY  DIFFUSION  COEFFICIENT  FOR  THE  MgCK-H-O  SYSTEM  (1  ATM) 
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AT  25  C  (1  ATM) 


curve  has  been  drawn  through  the  CaC^-I^O  data.  The  following  is  apparent: 

(i)  For  the  NaCl-l^O  system,  sufficient  data  is  available  at  all 

concentrations  only  for  25°C;  the  data  extrapolates  to  the  Nernst 
limiting  value;  the  Onsager-Fuoss  (52,  pp.  105-122,  243-255) 
theory  is  only  useful  for  the  dilute  region;  and  Gordon's  semi- 
theoretical  method  for  predicting  binary  ordinary  diffusion 
coefficients  in  concentrated  salt  solutions  yields  results  which 
are  quite  close  to  those  obtained  experimentally.  Gordon's 
method  requires  system  viscosity  and  partial  molal  volume  data 
for  the  temperatures  and  concentrations  of  interest.  Because  such 
information  is  generally  not  available,  this  method  is  of  limited 
utility. 

(ii)  For  the  MgC^-H^O  system,  very  little  data  at  any  temperature  or 
concentration  is  available;  no  effort  has  been  made  to  apply 
Gordon's  method  to  this  system;  and  the  Onsager-Fuoss  method 
appears  useful  in  the  very  dilute  region. 

(iii)  For  the  CaCl^-H^O  system,  only  data  at  25°C  is  available;  the 

unusual  variation  of  0^  with  concentration  at  this  temperature 
cannot  be  accounted  for  by  the  Onsager-Fuoss  theory;  if  a  vis¬ 
cosity  correction  is  made,  and  if  a  correction  is  made  for  the 
hydration  of  the  dissolved  salt,  then  the  shape  of  the  curve  at 
higher  concentrations  (including  the  maximum)  can  be  predicted  in 
an  approximate  manner  (66). 

In  Table  2.1  the  various  major  salt  constituents  of  a  typical  sea 
water  are  listed  along  with  the  associated  measured  values  of  and  o^ 
for  25°C  and  1  atm.  This  table  gives  an  indication  of  the  relative  orders 
of  magnitude  of  these  properties,  as  well  as  the  lack  of  essential  binary 
system  data.  Very  little  diffusion  data  has  appeared  in  the  literature 


*  The  data  given  by  Oholm  may  not  be  reliable. 
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RELATIVE  MAGNITUDES  OF  BINARY  SORET  COEFFICIENTS  AND  ORDINARY 
DIFFUSION  COEFFICIENTS  FOR  THE  MAJOR  SALT  CONSTITUENTS  OF 

SEA  WATER  (25°  C,  1  atm) 
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for  the  MgSO^-l^O  system,  while  none  has  been  published  for  the  CaSO^-I^O 
system  at  any  temperature.  Table  2.1  also  shows  the  relative  concentrations 
of  the  various  salt  constituents  of  ordinary  sea  water.  Most  of  the  binary 
and  o^2  are  reported  only  for  one  concentration,  that  is  0.01  molal  (gm 
moles  salt/1000  gm  water).  It  is  apparent  that  NaCl  is  by  far  the  major 
salt  constituent,  and  that  of  the  binary  diffusion  properties  which  are 
known  for  the  various  salts  in  sea  water,  the  (and  a^)  of  MgSO^-H^O 
and  the  D.„  of  KoS0, -H«0  are  much  different  than  the  D,„  for  the  other 
salt  systems  (at  0.01  molal).  No  data  indicating  the  effect  of  pressure  on 
D^2  for  any  of  the  binary  salt  solutions  listed  in  Table  2.1  has  been  re¬ 
ported. 

2.2  EXPERIMENTAL  PROGRAM 
2.2a  Experimental  Approach 

The  diaphragm  cell  method  was  selected  for  measuring  the  ordinary 
diffusion  coefficients  because  of  its  simplicity  and  accuracy  (approximately 
+  57.,  or  better)  (56,  pp.  253-261)  (49,  pp.  75-79). 

This  technique  was  developed  originally  by  Northrop  and  Anson  (91). 
Improvements  in  the  technique  have  been  made  by  Gordon  (92),  Stokes  (93) 

(60) (94) (95)  and  Robinson  (100).  Measurements  at  temperatures  other  than 
25°C  have  been  made  by  Firth  and  Tyrrell  (96).  Certain  modifications  of 
design  have  been  proposed,  notably  those  of  Lewis  (97)  (use  of  electrodes 
to  measure  differential  diffusion  coefficients  directly),  and  Smith  and 
Storrow  (98)  (use  of  a  Selsyn-transmitter *  for  stirring).  Direct  sampling 
during  the  diffusion  experiment  is  usually  not  employed,  though  such  methods 
have  been  used  (99). 

Methods  for  reducing  the  data  obtained  from  such  an  experiment 
have  been  given  by  Barnes  (101),  Gordon  (92),  Stokes  (60)(94)(95) ,  Firth 
and  Tyrrell  (96),  Dullien  and  Shemilt  (102) (103),  Olander  and  coworkers 
(51) (104) (105),  and  Robinson  (100). 


-14- 


The  general  feature  of  this  method  is  that  the  solute-solvent 
system  is,  at  the  start  of  the  experiment,  separated  by  a  porous  glass 
diaphragm  from  pure  solvent  in  the  upper  portion  of  a  vertical,  two- 
compartment  cell.  Stirrers  are  situated  on  both  surfaces  of  the  diaphragm 
so  that  all  the  solution  in  both  the  top  and  bottom  chambers  is  kept  at 
uniform  concentrations,  C_,  C  .  Thus,  molecular  diffusion  takes  place 
only  in  the  pores  of  the  diaphragm.  Because  it  is  impossible  to  predict 
the  diffusion  path  length  and  cross-sectional  area  for  such  a  system,  it 
is  necessary  to  calibrate  the  cell  with  a  solution  whose  ordinary  diffusion 
coefficient  is  known  at  the  temperature  of  interest  and  preferably  in  a 
concentration  range  similar  to  that  under  investigation.  Reliable  results 
have  been  obtained  at  concentrations  greater  than  0.05  molal  (>0.3  weight 
percent  NaCl  in  1^0,  for  example).  The  ordinary  diffusion  coefficients 
for  the  KCl-H^O  system  at  several  temperatures  and  concentrations  have 
been  accurately  determined,  and  KC1  concentrations 
have  been  accurately  determined  by  absolute  methods  (electrolytic 
conductance  and  optical).  £or  these  reasons,  this  particular  system  is 
useful  for  calibrating  the  cell  (determining  the  cell  constant)  (94). 

2.2b  Equipment  and  Procedures 

The  cell  used  to  measure  ordinary  diffusion  coefficients  was  a 
modification  of  the  diaphragm  cell  originally  developed  and  described  by 
Northrop  and  Anson  (91)  and  improved  by  Stokes  (93).  This  cell  provides 
data  free  from  errors  due  to  mixing  at  the  interface,  streaming  of  the 
heavier  liquid  through  the  diaphragm,  and  due  to  stagnant  layers  at  the 
surface  of  the  diaphragm.  In  addition,  our  modifications  of  the  cell 
permit  measurement  of  temperature  and  concentration  within  the  cell  at 
any  time  during  the  experiment. 

In  order  to  make  ordinary  diffusion  coefficient  measurements  at 
pressures  up  to  1000  bars  (~14,700  psi)  without  using  radioisotope  tracer 
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techniques,  a  new  cell  and  magnetic  stirring  system  had  to  be  designed 
which  would  have  operating  characteristics  similar  to  the  usual  glass 
atmospheric  pressure  diaphragm  diffusion  cells  and  yet  be  able  to  stand 
the  immense  pressures.  The  solution  developed  for  this  study  was  an 
ordinary  diffusion  bomb  of  the  design  shown  in  Figure  2.5.  The  fabricated 
bomb  is  shown  partially  disassembled  in  Figure  2.6.  K  Monel  was  chosen 
as  the  bomb  material  because  of  its  corrosion  resistance  to  saline 
solutions  (106),  its  nonmagnetic  property,  and  its  high  tensile  strength. 
For  the  pressure  range  of  the  experiments,  it  was  decided  that  unsupported 
area  seals  would  be  best  from  the  viewpoint  of  safety,  dependability,  and 
ease  of  assembly  and  disassembly  (107).  The  wall  thickness  of  the  bomb 
was  determined  by  using  an  average  of  estimates  calculated  from  three  con¬ 
servative  design  relations:  the  Maximum  Shear  Equation  (108),  the  ASME 
Code  Formula  (109),  and  the  Lame7  Formula  (109). 

One  unique  feature  of  this  diffusion  bomb  design  was  the 
inclusion  of  platinum  electrodes  attached  to  each  of  the  end  closure 
pieces  (see  Figures  2.6  and  2.8)  and  positioned  in  the  stirred  chambers 
on  either  side  of  the  porous  glass  diaphragm.  These  electrodes,  after 
proper  calibration,  could  be  used  to  follow  the  concentration  changes  in 
each  stirred  chamber  during  the  course  of  the  diffusion  process.  It  was 
originally  planned  that  by  obtaining  such  data  it  would  be  possible  to 
obtain  the  true  differential  diffusion  coefficient  directly  from  one 
experiment  without  having  to  make  a  number  of  experiments  at  several 
different  concentrations  and  then  treat  the  data  in  a  somewhat  involved 
manner  in  order  to  derive  the  differential  coefficients  from  the  measured 
integral  coefficients  (60)(94) (95) .  Originally  an  epoxy-type  unsupported 
area  seal  was  chosen  for  sealing  and  insulating  the  platinum  electrodes 
in  the  closure  pieces  because  of  its  apparent  simplicity,  high  electrical 


These  techniques  were  not  practical  for  the  salt  solutions  being 
investigated. 
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FIGURE  2.5.  ORDINARY  DIFFUSION  BOMB 


P03599 


U 


H 


sO 

CM 


o 

H 

Pm 


resistance,  and  corrosion  resistance  (110).  Problems  in  obtaining  a 
good  seal  between  deaerated  epoxy  cement  and  the  K  Monel  surface 
required  an  alternative  sealing  method.  The  one  finally  selected  was 
a  modification  of  the  usual  cone  -  pyrophyllite  (a  hydrous  aluminum 
silicate  -  Al^O^.  4Si02*  1^0)  electrical  lead  type  of  seal.  A  schematic 
representation  of  this  type  of  seal  is  shown  in  Figure  2.7. 

Magnetic  stirring  of  the  two  cell  chambers  on  either  side  of 
the  porous  glass  diaphragm  was  accomplished  by  using  epoxy  and  Teflon 
covered  magnets  supported  and  rotated  by  the  magnetic  stirring  system 
shown  in  Figure  2.8.  This  system  was  found  to  be  the  best  of  several 
designs  which  were  tried. 

Figure  2.8  also  shows  the  general  arrangement  of  the  bomb  in  the 
constant  temperature  bath  as  it  appeared  during  the  course  of  an  experiment. 
The  pressuring  system  shown  in  Figure  2.9  was  used  to  transmit,  measure, 
and  record  the  desired  pressure  in  the  bomb. 

It  was  found  that  the  most  reliable  conductance  measurements  of 

the  salt  solutions  in  the  two  bomb  chambers  were  made  when  there  was  no 

fluid  in  the  constant  temperature  bath  container.  This  was  because  there 

appeared  to  be  some  extraneous  and  erratic  capacitance  effects  between 

the  insulated  leads  to  the  platinum  electrodes  and  the  liquid  bath  medium 

in  contact  with  them.  A  Jones  bridge  (AC  resistance  bridge;  accuracy: 

4  7 

1  part  in  10  ;  sensitivity:  1  part  in  10  )  was  used  to  measure  the 
resistance  between  both  pairs  of  platinized  platinum  electrodes. 

Temperature  control  in  the  bath  (stirring  was  achieved  by  the 
rotating  magnetic  arms  shown  in  Figure  2.8)  was  maintained  to  within 
+  0.1°C  by  circulating  fluid  maintained  at  a  constant  temperature  (usually 
25.0°C)  through  the  helical  coil  situated  along  the  inside  wall  of  the 
bath  container. 

Salt  solution  concentrations  of  the  fluid  in  each  of  the  chambers 
were  measured  at  the  end  of  the  experiment  (after  depressurizing  to  1  atm, 
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SURFACES  EXPOSED 
TO  THE  SALT 


FIGURE  2.7.  MODIFIED  HIGH  PRESSURE  ELECTRICAL  SEAL  FOR 
ORDINARY  DIFFUSION  BOMB 
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FIGURE  2.8.  ARRANGEMENT  OF  ORDINARY  DIFFUSION  BOMB  IN  CONSTANT  TEMPERATURE  BATH 
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and  removal  of  the  solution  from  the  bomb)  by  the  use  of  a  small-volume 
(~2cc)  conductance  cell  (pipette-type  cell)  used  in  conjunction  with  the 
Jones  bridge  and  another  constant  temperature  bath. 

The  detailed  procedure  which  was  finally  developed  for  use  with 
this  bomb  is  given  in  Appendix  A.l 

The  volumes  of  both  the  top  and  bottom  chambers  and  the  diaphragm 
pores  were  measured  several  times  and  found  to  be  approximately  38.8  cc, 
39.1  cc,  and  0.5  cc,  respectively. 


2 . 3  RESULTS 


The  usual  measurements  made  in  an  ordinary  diffusion  cell 
experiment  are  the  molar  concentrations  of  the  solutions  in  each  chamber 
at  the  end  of  the  experiment,  the  temperature,  the  duration  of  the  experi¬ 
ment,  and  the  volumes  of  the  top  and  bottom  chambers  as  well  as  the 
diaphragm  pore  volume.  The  integral  binary  ordinary  diffusion  coefficient 
is  then  calculated  from  (51) (60) (92) (94-96)(100- 105) (assuming  no  volume 
change  on  mixing;  Vfi;  VDp//VT  as  0.1), 


D 


12 


Pt 


(12) 


where  the  cell  constant,  8,  is  defined  by 

*  ■  f  { V1  +  V1 }  (13> 

and  where  CT  is  usually  taken  as  0  in  experiments  with  binary  aqueous 
salt  solutions  since  pure  water  is  usually  put  into  the  top  chamber  at 
the  start  of  the  experiment. 

The  integral  coefficient  is  related  to  the  "true"  differential 
coefficient  (such  as  reported  in  Figures  2.2,  2.3,  and  2.4)  by  an  equation 


-23- 


of  a  form  such  that  x  in  concentration  regions  where  is 
independent  of  concentration.  This  is  the  case  for  NaCl-I^O  in  the 
vicinity  of  3.5  weight  percent  salt  at  25°C  and  1  atm. 

The  cell  constant  3  is  determined  by  performing  a  diffusion 
measurement  in  the  cell  with  a  solution  whose  is  well  known  (such 
as  KC1-H20  at  25°C). 

The  actual  C  is  not  the  molar  concentration  of  the  salt 

"o  ‘ 

solution  originally  put  into  the  bottom  chamber  but,  because  a  pre¬ 
liminary  diffusion  experiment  is  always  allowed  to  take  place  before  the 
actual  experiment  is  started  (in  order  to  establish  an  approximately 
linear  concentration  gradient  through  the  diaphragm  before  the  start  of 
the  experiment),  it  must  be  calculated  by  a  mole  balance  from, 


[< 


VT  +  *VDP  I 
VB  +  *VDP  J 


(14) 


The  volume  ratio  was  nearly  equal  to  unity  (within  0.47.)  for  the  geometry 

employed  in  the  ordinary  diffusion  bomb. 

The  measurements  made  with  the  ordinary  diffusion  bomb  during 

£ 

the  first  phase  of  the  experiments  did  not  make  use  of  the  electrodes. 

The  salt  concentrations  were  obtained  after  the  bomb  had  been  brought 
back  to  1  atm  by  removing  samples  and  making  use  of  a  small  pipette  con- 
ductivity  cell  and  the  Jones  bridge.  The  concentrations  measured  in 
this  manner  are  reported  (as  C_  and  C_  )  along  with  other  pertinent 

ifi  Bf • 

associated  data  in  Table  2.2.  The  overall  duration  of  the  experiment, 
t',  and  the  "back-calculated"  initial  bottom  chamber  concentration  (C^,  =  o 


non 
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'B, 


+  c„ 


o 

(14a) 


*  Initially  some  problems  were  encountered  with  their  performance. 
**  Industrial  Instruments,  Inc. 


ORDINARY  DIFTIJS ION  BCHB  EXPERIMENTS  —  FIRST  PHASE 
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0.193  (cm  )  Dia.  #  III  (Experiments  0DB-3B  to  ODB-32B) 

-  2 

0.2146  (cm  )  Dia.  #  IV  (Experiments  ODB-33B  to  ODB-43B) 


were  used  to  calculate 


(  St')"1 


(12a) 


which  are  also  tabulated  in  Table  2.2.  The  3  values  were  determined  by 
making  diffusion  measurements  with  KCl-l^O  solutions.  The  considerable 
scatter  in  the  Dj2  data  obtained  for  experiments  run  under  seemingly 
identical  conditions  (compare  the  Dj2  data  from  Experiments  ODB-28B  to 
ODB-35B)  suggested  that  some  random  processes  or  effects  were  significantly 
altering  the  results  from  one  experiment  to  the  next  despite  the  care  which 
was  taken  to  duplicate  experimental  conditions  and  techniques. 

Based  on  these  results,  a  second  phase  of  experiments  was  per¬ 
formed  in  which  a  systematic  effort  was  made  to  uncover  the  exact  source 
of  the  data  scatter.  The  results  obtained  from  this  second  phase  of 
experiments  are  reported  in  Table  2.3  along  with  the  experimental  (and 
data  reduction)  conditions  peculiar  to  each  experiment. 

A  differential  error  analysis  of  the  effect  of  pressure  on  the 
diaphragm  cell  constant,  3»  revealed  that, 

A3  _  AA  A 1  ,  AV  nex 

3  "a  "  1  V  v  ' 

For  the  porous  glass  disc  used  (Kimax  borosilicate  porous  discs),  A 1 / 1 

was  estimated  to  be  less  than  0.01%  in  going  from  1  atm  to  1000  atm  (1 

AA 

taken  to  be  a  constant  fraction  of  the  diaphragm  thickness),  while 
was  estimated  to  be  less  than  0.2%  under  the  same  conditions.  The  quantity 
AV/V  is  best  estimated  by  direct  experimental  measurement  (112)  but  would 
be  expected  to  be  at  least  less  than  5%  for  these  conditions.  Thus  changes 
in  3  with  pressure  cannot  account  for  the  data  scatter. 
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tail*  a. i 

ordinary  Dirrmiow  km  nmnnri  -  ncow  nuu 


ixnRiMirr 

NIMIR 

0DIRS' 

•ALT  IN  AQUIO'J* 
SOLUTION 

DURATION  Of 
IXmilCNT 

(«0,-tf,) 

(aacaada) 

V 

[weight 
L  par  can 

V 
J  RSS 

V 

fwaight  1 
LpereentJ 

V 

[W& 

■V1# 

f weight  * 
i  percent. 

lldkll 

ir«is] 

DIAMRA0H  6 

RUM! 

(ca‘J) 

5;,  a  10* 

(ca*/aac) 

1 

Ml 

7,200 

0.184 

0.031 

3.340 

0.381 

3.324 

0.612 

VII 

-- 

*• 

2 

ruci 

7,200 

0.112 

0.018 

3.273 

0.570 

3.387 

0.366 

VII 

-- 

-• 

3 

NaCl 

7.200 

0.0)3 

0.005 

3.515 

0.613 

3.548 

0.618 

vu 

0.16)6 

1.2)0 

4 

KCt 

94,300 

0.438 

0.059 

3.049 

0.416 

3.487 

0.4  73 

VII 

0.1836 

5 

KC1 

71,520 

1  o  286 

0.17) 

3.263 

0.446 

4.549 

0.619 

VII 

0.6213 

6 

NaCl 

7,200 

0.133 

0  025 

2.920 

0.507 

3.075 

0.3)2 

VII 

-- 

-- 

7 

KC1 

57,7*4 

0.  108 

0.014 

3.625 

0.324 

3.433 

0.538 

VII 

0.05012 

8 

NaCl 

68,400 

0.622 

0.1063 

3.068 

0.333 

3.690 

0.63*5 

VII 

0.16)6 

1.256 

9 

NaCl 

73,320 

0.434 

0.078 

3.335 

0.581 

3.7*4 

0.639 

VII 

0. 18)6 

1.346 

to 

NaCl 

73,260 

0.518 

0.084 

3.520 

0.614 

4.0)8 

0.70) 

VII 

0. 1836 

1.343 

u 

XC1 

242,040 

1.800 

0.244 

2.193 

0.297 

).«*) 

0.541 

VII 

0.3220 

•• 

12 

NaCl 

2,400 

0.105 

0.018 

-• 

-* 

-- 

-• 

VIII 

- 

•• 

13 

NaCl 

54,784 

0.345 

0.058 

2.340 

0.441 

2.885 

0.4*9* 

VIII 

0.2525 

1.75 

14 

RC1 

2)4,876 

1.462 

0.1 175 

2.412 

0.347 

4.374 

0.5*45 

VIII 

0.2325 

•• 

COMNTI 


Pmiurmt  Ion  arror  dataralnat  Ion;  IJ 
rlnutaa  to  praaaurlaa  to  1000  kara,  U  houro 
•t  1000  kora,  15  alnutaa  to  dapraaaurlta 

Praaaur Itat Ion  arror  dataralnat Ion;  1  hour 
lo  praaaurlaa  to  1000  kara,  1  hour  to 
dapraaaurlta 

Oalarolnat  I  on  of  tha  axtant  of  ordinary 
dlffualon  at  1  bar  for  a  2  hour  parlod 

Calibration  of  Dlaphraga  VII  for  5 

Calibration  of  Dlaphrai*  VII  for  8;anail  crack 
In  Dlaphrafo  probably  praaant  bafora  aaparl—a—t 
Praaaur I aat ion  arror  dataralnat Ion;  15 
■  lnutaa  to  praaaurlaa  to  1000  kora,  lb  houra 
at  1000  kara.  15  alnutaa  to  dapraaaurlta 

Calibration  of  Dlaphraga  Vll  for  6 

Ordinary  dlffualon  aaparlMnt  at  1000  kara; 

15  alnutaa  to  praaaurlaa,  15  alnutaa  to 
dapraaaurlta 

Ordinary  dlffualon  axparlaant  at  1000  bora; 

15  alnutaa  to  praaaurlaa,  15  alnutaa  to 
dapraaaurlta 

Ordinary  dlffualon  aapar<aant  at  1000  kara; 

15  alnutaa  to  praaaurlaa,  15  alnutaa  to 
dapraaaurlta 

Calibration  of  Diaphrafa  VII  for  8;  flna 
crack  la  dlaphraga  dlacovarad  at  an!  of 
axparlaant 

Vaa  of  top  alactroda  pair  in  boa*  to  chock 
praaaurltat Ion  arror;  boa*  la  air  bath;  no 
at Irrlng  for  2  houra  Initially,  than 
praaaur laad  to  1000  bara  In  ~15  -lnutaa, 
bald  at  noo  bara  for  -10  alautaa,  and  than 
Rapraaaiiriaad  la -15  alautaa;  la  a  tu 
coacaatrat Iona  aaaaurad  vltb  tap  alaetrodaa 

Ordinary  dlffualon  aaparlaant  at  1000  kara; 


top  alaetrodaa  uaad,  concantrat Iona  at 
datarnlaad  ulth  plpatta  call 

Calibration  of  Olaphrapha  VIII  for  ( 


‘f 


tXPIklWlfl  SALT  IH  AQUEOUS 
NTMIR  SOLUTION 
O0HS- 


U* 


NaCl 


13** 


NaCl 


15 


16 


NaCl 


NaCl 


N0TI:  T  -  25°C 


DURATION  OP 
Einintn 

<V*f> 

(eeconda) 

57,480 


57,490 


V 


V. 


feiSJ!]  fSlU|  fsaissl  fSEliS  i  rssilil  faclat!  r-olaal  r-olaal 

LlltarJ  t  1  ltarj  LlitarJ  LlitarJ  1  UtarJ  tlltarj  LTh^j  iTuTJJ 


PUntRACN 

NIMIR 


0,2  a  105 


0.058 


0.058 


0.550 


0.061 


0.46 1 


0.031 


0.464 


0.008 


0.4*2 


0.441 


0.0663  0.4*6 


0.0183 


0.544 


0.330 


0.061 


0.3*4 


0.0)3 


0.3*4 


0.0)3 


(«  l)  (cn^/aae) 

VIII  0.2525  0.144  Ixparlaant  008RS-I3  vlth 

actual  tlaa  (t  -t  )  and  con- 
cant rat  Iona  uadd  to  calculata 
"  flfit  calibration  curve 

®*  *ta*  alaetrodaa  v»-  c  “"d 
(Cg,  uiad  for  Cl;  Cj  ,C( 

datarnlnad  by  sola  bilanffaa 

VIII  0.2525  1.3*  Ixparlaant  00RRI-13  with 

actual  tlaa  <t#-tf>  and  con- 
£ant rations  uaaa  to  calculata 
Oil  —  aacaad  calibration 

‘/"L",**'9’  V*- 

C  uaad  (Cf.  uaad  for  C>; 

^  iataralnad  by  nolo 

bafaacas 

Calibration  of  top  alaetrodaa 
1*  bo*;  concantrat  Iona 
dataratnod  by  plpatta  call; 

P  •  1000  bara 

Calibration  of  top  alaetrodaa 
In  book;  concantrat Iona 
datarnlnad  by  plpotta  call; 
t  •  1000  bara 
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Since  the  molar  concentrations  that  were  measured  at  1  atm 
pressure  in  the  pipette  cell  enter  into  the  calculation  of  38  3 
ratio  (see  Eqs.  12  or  12a),  the  effect  of  pressure  on  the  C's  can  be 
neglected  over  the  pressure  range  of  interest. 

An  analysis  of  random  errors  in  the  measurement  of  t  and  C's 
and  their  effect  on  shows  that  (93), 


AD 


12 


D 


12 


(16) 


where  f  is  a  function  of  C  and  C  and  increases  as  t,  increases.  This 

.  f  ®f  l 

function  always  had  a  value  of  less  than  1.5  for  the  experiments  reported 
here.  The  experiment  duration  could  be  determined  to  better  than  0.1% 
while  the  accuracy  of  the  concentration  measurements  using  the  extensively 
calibrated  pipette  conductivity  cell  was  better  than  two  percent. 

These  possible  sources  of  error  could  not,  therefore,  account 
for  the  magnitude  of  the  data  scatter  reported  in  Table  2.2. 

The  first  results  reported  in  Table  2.3  reveal  that  if  the 
ordinary  diffusion  bomb  is  brought  up  to  pressure  and  depressurized 
very  slowly,  less  nondiffusive  bulk  transport  occurs  across  the  diaphragm. 
Repeated  runs  under  nearly  identical  conditions  could  give  nearly  the  same 
results  (Experiments  0DBRS-9  and  10).  The  porous  glass  diaphragm  was  sus¬ 
ceptible  to  cracking  as  was  found  after  Experiment  0DBRS-11. 

The  evidence  thus  obtained  suggested  that  a  nondiffusive  bulk 
transport  of  solution  occurred  during  both  the  pressurizing  and  depressurizing 
processes  even  though  pressure  was  simultaneously  transmitted  to  both  the 
top  and  bottom  chambers  of  the  bomb  through  high  pressure  stainless  steel 
capillary  tubing  of  identical  lengths  (pressure  drops  equal). 
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The  use  of  the  electrodes  attached  to  the  end  closure  pieces 
was  then  initiated.  Early  problems  encountered  with  these  electrodes 
were  at  least  partially  overcome  by  not  using  fluid  in  the  constant  temp¬ 
erature  bath  (a  crude  constant  temperature  air  bath  was  used  instead), 
by  the  use  of  an  improved  method  for  attaching  the  lead  wires  to  the 
platinum  electrode  wires,  and  by  improving  the  operating  techniques  used 
to  obtain  transient  resistance  values  with  the  Jones  bridge. 

The  calibration  of  the  pair  of  electrodes  in  the  top  chamber  (those 
in  the  bottom  chamber  were  inoperative  due  to  a  broken  lead  wire)  is  shown 
in  Figure  2.10  for  the  pressure  limits  of  interest.  The  data  shown  graph¬ 
ically  in  Figures  2.11  and  2.12  vividly  show  the  presence  of  the  bulk 
transport  effect  on  pressurizing  and  depressurizing  as  well  as  its  magnitude 
and  randomness.  In  Figure  2.12,  the  majority  of  the  ordinary  diffusion 
process  occurred  between  tQ  and  tf.  The  pressurizing  time,  tQ,  to  tQ,  and 
the  depressurizing  time,  t^  to  t^ ,  ,  were  kept  small.  The  concentrations  at 
t_,  were  obtained  using  the  pipette  cell.  C  and  C  were  measured  in 
situ  with  the  calibrated  electrodes.  C_  andrC  wer8  calculated  according 

Of  B 

to  the  scheme  outlined  in  Appendix  A. 2.  r0DBRS-l§*  denotes  the  case  when 

the  measured  top  electrode  pair  resistances  from  Experiment  ODBRS-13 

were  reduced  to  molar  concentrations  by  the  first  calibration  (see  Figure  2.10), 

while  ODBRS**  denotes  the  case  when  the  second  (more  valid)  calibration  was 

used. 

A  summary  of  essentially  all  the  experimentally  measured  and 
data  are  shown  in  Figure  2.13  for  the  various  pressures  studied.  As 
already  noted,  the  scatter  of  the  data  is  considerable.  If  the  data 
from  ODBRS-13  is  treated  in  the  same  way  as  was  the  data  obtained  from 
the  other  experiments,  the  shown  above  the  mid-scatter  line  would  be 
obtained.  However  if  the  electrode  measurements  from  this  same  experiment 
were  employed  to  determine  in  situ  concentration  measurements  at  all  times 
during  the  experiment,  then  the  interesting  result  shown  in  Figure  2.13 


-29- 


10/  /R  (p,  mho) 


C  (MOLES/LITER) 

FIGURE  2.10.  TOP  ELECTRODE  PAIR  CALIBRATION 


PI 1042 
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PI 1039 

FIGURE  2.11.  VARIATION  OF  CONDUCTANCE  WITH  PRESSURE  IN  TOP  CHAMBER  OF 

ORDINARY  DIFFUSION  BOMB 
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(MOLES /LITER) 


P11043 


FIGURE  2.12.  VARIATION  OF  MOLAR  CONCENTRATION  WITH  TIME  IN  THE 

ORDINARY  DIFFUSION  BOMB 
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Integral  binary  diffusion  coefficient 

F  x  I05  (CM2  /SEC) 


2.5 


DATA  [ONLY  C.  C_  MEASURED;  C„  E  ol 

L  V'  Tf'  T0  J 


DATA  -  EXPERIMENTAL  LIMITATIONS  foNLY  C.  C-  MEASURED;  C„  E  ol 

L  V*  V  T0  J 

LITERATURE  VALUE 

EXPERIMENT  ODBRS -13**  TlN  SITU  MEASUREMENT  OF  CONCENTRATION  IN  TOp"| 

[CHAMBER  ONLY 


J  EXPERIMENT  ODBRS-13**  [ IN  SITU  MEASUREMENT  OF  CONCENTRATION  IN  TOpI 

[chamber  only  J 

)  VALUE  BASED  ON  EXTRAPOLATION  OF  A  (?)  DATA  USING  A H  •  CONSTANT 

(WALDEN'S  RULE) 

}  VALUE  BASED  ON  EXTRAPOLATION  OF  A(?)  DATA  USING  AM?-CONSTANT 

(MODIFIED  WALDEN'S  RULE) 

NUMBERS  DENOTE  EXPERIMENT  (SUFFIX  B  OMITTED)  (ODB  EXPERIMENTS  UNLESS 

OTHERWISE  INDICATED) 

DOTS  DENOTE  EXPERIMENTS  OF  LONGER  TIME  DURATION  (  >  40  HOURS)  -  OTHER 

DATA  ARE  FOR  EXPERIMENTS  OF  <  28  HOURS  DURATION 

FLAGS  DENOTE  EXPERIMENTS  PERFORMED  WITH  INCREASED  STIRRING  RATE  (70  RPM) 

- DENOTES  APPROXIMATE  LIMITS  OF  SCATT.A  OF  THE  MAJORITY  OF  DATA 

■  —  DENOTES  APPROXIMATE  AVERAGE  OF  DATA  SCATTER  LIMITS 


s*' 

i»n 

32  CT 


ODBRS 


-13 


012  / 

©14,17 


34^ 


©22 

"A 

25® 

ODBRS -13** 

06 _ 

*5Cf 

400  500  600 

PRESSURE  (BARS) 


PI 1030 


FIGURE  2.13.  INTEGRAL  BINARY  DIFFUSION  COEFFICIENT  OF  3.5  WEIGHT 

PERCENT  AQUEOUS  SODIUM  CHLORIDE  SOLUTION  AT  25°C 
AS  A  FUNCTION  OF  PRESSURE 
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(large  lowest  hexagon)  is  found  (second  calibration  used  —  see  Appendix 
A. 2). 

Assuming  the  validity  of  the  Stokes-Einstein  expression  (see 

Eq.  9)  as  a  first  approximation,  a  crude  estimate  may  be  made  of  the 

effect  of  pressure  on  (°r  D^)  using  available  conductance  data 

obtained  at  high  pressure  ( 1 13 ) ( 1 14 )  and  either  Walden's  rule  (p,  A  = 

constant)  or  a  modified  form  of  Walden's  rule  which  appears  to  be  more 

2 

valid  for  salt  solutions  (p,  A  =  constant)  (113).  The  results  of  estimates 
of  Dj2  at  1000  bars  based  on  both  of  these  extrapolation  methods  are  shown 
in  Figure  2.13. 

From  these  few  results  it  may  be  now  concluded  that  pressure 
seems  to  have  little  effect  on  the  for  nominally  3.5  weight  percent 
NaCl  in  1^0  at  25°C  over  the  range  1-1000  bars.  Whether  the  effect  is 
to  slightly  increase  or  slightly  decrease  at  1000  bars  compared  with 
its  value  at  1  bar  (for  25°C)  remains  to  be  determined  by  further  experi¬ 
mentation. 

Future  ordinary  diffusion  bomb  experiments  (with  binary  solutions) 
conducted  at  pressures  greater  than  1  bar  must  use  both  top  and  both 
electrode  pairs  which  have  been  extensively  calibrated.  Further,  care  must 
be  exercised  that  the  pressure  transmitting  lines  to  the  bomb  are  flushed 
clean  before  each  experiment.  Under  these  conditions,  can  be  cal¬ 
culated  with  confidence  directly  from  Equation  12. 

This  new  experimental  technique  (use  of  the  ordinary  diffusion 
bomb)  appears  to  offer  a  number  of  advantages  over  radioisotope  techniques 
(115-118),  especially  for  binary  salt  systems  where  in  situ  concentration 
measurements  are  possible  without  incurring  radiation  hazards. 

2.4  EXTENSION  TO  MULTICOMPONENT  SYSTEMS 

The  theoretical  description  of  ion  diffusion  under  isothermal 
conditions  in  a  multicomponent  system  has  been  presented  in  the  literature 
(55)  (119-121).  A  basic  result  is  that,  to  a  good  approximation,  multi- 
component  ionic  ordinary  diffusion  coefficients  can  be  measured  with  the 
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diaphragm  diffusion  cell  using  Eq .  12,  but  interpreting  as  Di  (a 

multicomponent  ionic  ordinary  diffusion  coefficient  dependent,  in 
general,  upon  temperature,  pressure,  and  the  relative  concentrations  of  the 
ions  present)  and  interpreting  the  molar  concentrations  as  ionic  molar 
concentrations  of  the  species  i. 

Since  conductance  measurements  in  multicomponent  salt  solutions 
do  not  yield  concentration  values  for  the  individual  ions,  some  other  in 
situ  concentration  method  must  be  sought  (such  as  spectrophotometry  using 
high  pressure  optical  windows  in  the  bomb  walls)  or  displacement  sampling 
under  pressure.  For  the  types  of  salts  of  most  abundance  in  sea  water, 
the  high  pressure  sampling  technique  would  probably  be  the  most  practical. 
Mercury  could  be  used  as  the  displacing  fluid.  High  pressure  valves 
could  be  attached  to  the  side  wall  of  the  ordinary  diffusion  bomb  at 
points  both  above  and  below  the  diaphragm.  The  downstream  side  of  these 
valves  could  be  pressurized  to  a  pressure  just  slightly  greater  than  the 
bomb  pressure  before  sampling  would  be  initiated.  With  the  stirrers 
stopped,  the  valves  could  be  opened  and  then  the  mercury  could  be  injected 
at  the  same  time  that  a  sample  is  withdrawn  into  the  tubing  attached  to 
the  downstream  side  of  the  valves  with  the  pressure  on  this  side  of  the 
valves  being  maintained  constant  (by  increasing  the  volume  of  the  system, 
i.e.,  by  creating  sample  volumes  in  the  lines).  After  having  "trapped" 
samples  at  high  pressure  under  conditions  such  that  bulk  flow  across  the 
diaphragm  would  be  minimum,  the  pressure  on  the  contained  and  isolated 
samples  could  be  reduced  to  1  bar.  Individual  ion  concentrations  could 
then  be  measured  by  flame  photometric  techniques. 
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SECTION  3 


SORET  COEFFICIENT  STUDY 


3.1  ANALYTICAL  CONSIDERATIONS 

The  phenomenon  of  thermal  diffusion  in  salt  solutions  has  been 
recently  reviewed  by  Agar  (70)  and  Tyrrell  (49).  There  hay  been  no 
attempt  to  review  the  binary  Soret  coefficient  data  for  many  of  the  salts 
in  sea  water,*  nor  to  apply  considerations  of  multicomponent  effects  to 
synthetic  or  actual  seawater  (55). 

The  Soret  coefficient  of  dilute  and  concentrated  binary  salt 
solutions  has  been  of  interest  for  some  time.  Early  work  by  Tanner  (72) 

(73) (74) (75)  and  more  recent  work  by  Agar  and  coworkers  (76) (70) (77) (78) 
(79)(80)(81)(82)(61,  pp.  200-223),  Turner  and  coworkers  (68)(69),  Chanu 
and  coworkers  (83)  (84) (85) (86) ,  Longsworth  (87)  and  Thomaes  (88)  has 
consisted  primarily  of  experimental  measurements  of  the  binary  Soret 
coefficient,  optical  or  conductivity  methods. 

The  only  theoretical  attempts  which  have  been  made  in  order 
to  predict  binary  Soret  coefficients  of  salt  solutions  have  been  based 

Von  Halle  (71)  has  given  a  tabular  review  of  all  binary  Soret  coefficient 
data  published  before  1957. 

**  The  theory  of  Bearman  (89)  is  not  applicable  to  salt  solutions. 
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on  molar  heat  of  transfer  data  for  the  individual  ions  at  the  temperatures 
and  concentrations  of  interest  (49).  These  ionic  molar  heats  of  transfer 
can  be  related  to  corresponding  quantities  for  the  infinitely  dilute  state. 

No  theoretical  means  are  available  for  calculating  these  latter  quantities. 

Ionic  molar  heats  of  transfer  are  usually  determined  indirectly  by  measuring 
the  binary  Soret  coefficient.  Thus,  in  effect,  there  are  no  theoretical 
tools  by  which  cr^  values  may  be  calculated. 

The  major  salt  constituents  of  a  typical  sea  water  and  their 
associated  Soret  coefficients  for  one  concentration  (at  25°C)  are  listed 
in  Table  2.1.  From  this  table  of  values,  it  is  apparent  that  significant 
differences  exist  between  the  various  values  of  data  at  any 

temperature  or  composition  is  available  for  those  salts  for  which  o^  values 
are  not  listed.  Even  in  the  case  of  those  for  which  a  ^  values  are  given, 
little  or  no  reliable  data  exists  in  the  literature  for  the  temperatures 
and  concentrations  of  interest  here.  No  data  indicating  the  effect  of 
pressure  in  cr^  f°r  aqueous  salt  solutions  has  been  published. 

All  published,  reliable,  Soret  coefficient  data  for  the  NaCl-I^O 
and  MgC^-H^O  systems  as  a  function  of  composition  and  temperature  are  pre- 
sented  in  Figures  3.1  and  3.2.  Smooth  curves  have  been  drawn  through  some 
of  the  data.  The  following  is  apparent: 

(i)  For  the  NaCl-I^O  system,  data  at  temperatures  other  than  25°C 
is  lacking;  data  in  the  3.5  weight  percent  region  even  at  25°C 
is  not  available;  recent  data  obtained  by  Price  (79)  suggests 
that  0^2  does  go  through  a  minimum  as  the  concentration  is  increased; 
the  only  data  available  at  a  temperature  other  than  25°C  indicates 
that  the  temperature  dependence  of  the  Soret  coefficient  is  appreciable. 

* 

The  sources  for  the  data  are  as  follows:  a-  (NaCl-H«0):  Agar  (61),  pp. 
200-223),  Agar  and  Turner  (78),  Snowdon  and  Turner  (68) ,  Price  (79),  Chanu 
(83)  (85)  and,  Chanu  and  Lenoble  (83) 

<7^2  (MgC^-l^O) :  Payton  and  Turner  (69),  Cova  (90),  and  Tanner  (75). 
a  12  (KCI-H2O):  Longsworth  (106). 
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FIGURE  3.2.  SORET  COEFFICIENTS  FOR  THE  MgClj-I^O  SYSTEM  (1  ATM) 


(ii)  For  the  MgCl-i^O  system,  essentially  no  data  at  any  temperature 
is  available  over  a  large  concentration  range;  the  Soret 
coefficient  is  slightly  negat ive  at  high  concentrations  (salt 
migrates  to  the  hot  wall);  the  Soret  coefficient,  over  most 
of  the  concentration  range,  is  significantly  different  than 
that  of  the  NaCl-l^O  system  at  the  same  temperature. 

It  is  of  interest  to  consider  the  temperature  dependence  of  the 
binary  Soret  coefficient.  This  is  shown  in  Figures  3.1,  3.3,  and  3.4  for 
the  NaCl-I^O  and  KCl-I^O  systems.  Longsworth's  data  are  almost  the  only 
data  which  have  been  reported  that  consider,  systematically,  the  temperature 
dependence  of  a ^  for  a  binary  salt  solution  (in  this  case,  KCl-H^O). 

Figure  3.4  clearly  shows  the  strong  temperature  dependence  of  o^  for  both 
dilute  and  concentrated  solutions. 

3.2  EXPERIMENTAL  PROGRAM 
3.2a  Experimental  Approach 

After  surveying  the  various  methods  available  for  measuring  the 
Soret  coefficient  of  binary  salt  solutions,  it  was  apparent  that  a  con¬ 
duct  ime  trie  method  offered  the  only  single  technique  capable  of  making 
measurements  in  both  dilute  and  concentrated  solutions  (optical  methods 
are  limited  to  relatively  concentrated  solutions  in  which  index  of  re¬ 
fraction  gradients  are  large  enough  to  be  readily  measured)(87) .  Of  the 
two  principal  conductimetz ic  methods  (77) (78),  the  cell  design  developed 
by  Agar  and  Turner  (77)  seemed  to  offer  the  advantages  of  simplicity,  and 
ease  in  reduction  of  the  data.  In  addition,  there  appeared  to  be  fewer 
assumptions  and  uncertainties  associated  with  the  theory  of  the  cell 
operation.  Agar  and  Turner  (77)  had  suggested  that  this  cell  would  be 
useful  at  higher  concentrations  (>0.05  molal)  where  convective  mixing 
would  be  expected  to  be  less  of  a  problem  (the  hot  wall  is  always  situated 

on  top).  This  technique  required  the  use  of  an  AC  electrolytic  conductivity 

4 

bridge  sensitive  to  at  least  1  part  in  10  ,  and  required  the  temperature  at 
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FIGURE  i.4  COMPARISON  OF  THE  TEMPERATURE  DEPENDENCE  OF  THE  BINARY  SORET 
COEFFICIENTS  OF  THE  NaCl-H-0  and  KCl-Ho0  SYSTEMS  (1  ATM) 


the  hot  and  cold  walls  of  the  cell  to  remain  constant  within  0.01°C 
throughout  the  entire  duration  of  the  run  (runs  lasted  1  to  2  days). 

3.2b  Equipment  and  Procedures 

The  Soret  cell  used  during  these  tests  was  a  modification  of 
that  designed  and  used  by  Agar  and  Turner  (77).  The  cell  basically 
consists  of  two  jackets  for  heating  and  cooling  with  a  chamber  to  contain 
the  solution  between  them.  Provision  is  made  to  measure  temperatures  and 
concentrations  as  a  function  of  time. 

The  unique  feature  of  the  cell  used  in  this  experimental  program 
is  the  fact  that  it  can  be  used  to  make  thermal  diffusion  measurements  at 
pressures  up  to  1000  bars.  A  simplified  schematic  illustration  of  the 
cell,  called  a  Soret  bomb,  is  given  in  Figure  3.5.  The  detailed  design 
of  the  bomb  and  the  associated  heating  and  cooling  chambers  (through 
which  constant  temperature  fluid  is  circulated)  is  shown  in  Figure  3.6.  The 
size  of  the  central  chamber  and  the  mode  of  installation  of  the  electrode 
seals  are  shown  in  Figure  3.7.  As  in  the  case  of  the  ordinary  diffusion 
bomb,  the  central  cell  section  and  the  end  plates  of  the  Soret  bomb  were 
fabricated  from  K  Monel.  The  wall  thicknesses  were  determined  from  the 
usual  pressure  vessel  design  equations  (108)  (109).  The  bolts  which  hold 
the  two  end  plates  in  position  also  are  insulated  from  the  rest  of  the 
bomb  in  order  to  minimize  nonuniform  heat  transfer  problems  in  the  central 
cell  body  itself.  The  size  and  number  of  bolts  necessary  for  the  pressures 
of  interest  were  determined  using  known  techniques  of  fastener  design  (111). 

Whereas  unsupported  area  seals  were  used  in  the  ordinary  diffusion 
bomb,  in  the  Soret  bomb  o-ring  seals  with  back  up  rings  were  used.  This 
difference  is  due  to  the  requirement  that  the  heat  flux  through  the  central 
hole  region  of  the  Soret  bomb  wherein  the  salt  solution  is  contained  should 
be  uniform  and  one-dimensional.  It  is  therefore  undesirable  to  use  bulky 
gasketing  materials  having  thermal  conductivities  much  different  than  the 
surrounding  medium.  Phenolic  resin  plates  were  used  in  the  Soret  bomb 
in  order  to  make  the  thermal  resistances  (to  heat  transfer)  nearly  equal 
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FIGURE  3.5.  SORET  BOMB  ARRANGEMENT  (CELL  FASTENING  BOLTS  NOT  SHOWN  -  ENTIRE  ASSEMBLY 
ENCLOSED  IN  A  CONSTANT  TEMPERATURE  AIR  BATH) 
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FIGURE  3.7.  SORET  BOMB  DETAILS  —  CONDUCTIVITY  PROBES  AND  PRESSURE 


through  each  section.  Referring  to  Figure  3.8,  it  is  evident  that  if  the 
thermal  resistance  along  path  A-A*  is  to  be  the  same  as  along  B-B',  then 
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that  is,  the  thickness  of  the  phenolic  resin  plates  could  be  determined 
after  the  cell  height  had  been  specified. 

The  size  of  the  central  cavity  in  the  Soret  bomb  was  made  similar 
to  that  used  by  Agar  and  Turner  (77). 

The  operation  of  the  Soret  bomb  was  carried  out  in  a  manner  exactly 
similar  to  that  used  by  Agar  and  Turner  (77),  except  that  the  inside  of  the 
bomb  was  pressured  up  to  the  desired  pressure  before  the  start  of  an  experi¬ 
ment  by  the  transmitting  system  shown  in  Figure  2.9.  It  was  found  to  be 
especially  important  to  maintain  the  ".emoerature  of  the  heated  and  cooling 
circulating  fluids  (an  oil)  constant  to  within  0.0X°C  throughout  the 
duration  of  the  experiment.  In  addition,  very  careful  use  of  the  Jones 
bridge  was  required  (input  bridge  voltage  ~  0.25  volts;  voltage  frequency 
~  4000  cycles/sec)  in  order  that  the  appropriate  accuracy  in  resistance 
measurements  could  be  obtained.  Resistance  changes  during  the  course  of 
the  experiment  were  not  greater  than  about  2  ohms. 


3.3  RESULTS 

By  carefully  controlling  the  temperatures  of  the  heated  upper 
p^-te  and  the  cooled  lower  plate  to  within  0.01°C  over  a  twenty  four  hour 
period,  and  by  exercising  caution  in  the  use  of  the  highly  accurate  Jones 
bridge,  it  is  possible  to  make  reliable  transient  resistance  measurements 
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(with  the  electrodes  in  the  Soret  bomb)  which  are  indicative  of  the 
thermal  and  ordinary  diffusion  processes  occurring  in  the  center  cavity 
of  the  Soret  bomb.  The  binary  Soret  coefficient  is  then  calculated  from 
(for  t  >  0/3)  (111), 


R  -  R 
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r  x  i  ■ 
o  L  a  '  2  _ 


(18) 


where  AT  represents  the  temperature  difference  between  the  inside  surfaces 
of  the  top  and  bottom  end  plates  and  R  and  R  are  the  stationary  state 
and  initial  electrode  pair  resistances  determined  by  appropriate  extra¬ 
polation  of  the  resistance- time  data. 

After  making  a  considerable  number  of  improvements  in  the 
experimental  apparatus  (baths,  bomb,  and  bridge  connections)  it  was  finally 
possible  to  make  the  first  reliable  measurement  of  o^  for  the  NaCl-l^O 
system  at  a  concentration  similar  to  that  of  sea  water  and  at  a  pressure 
of  approximately  1000  bars  (T~27°C)  (Experiment  SDB-6:  top  electrodes). 

The  raw  and  derived  data  are  tabulated  in  Table  3.1.  Figures 
3.9  and  3.10  show  how  the  data  was  treated  in  order  to  obtain  the  required 
extrapolated  quantities,  R  and  R  .  The  consistency  and  reliability  of 

CD  O 

the  data  is  indicated  in  these  figures  (large  time  data  (t  >  0/3)  is 
weighted  most  strongly).  The  conformity  of  the  data  to  straight  lines, 
as  shown  in  these  figures,  is  required  for  consistency  with  the  assumptions 
implicit  in  Eq.  18  (111). 

Figure  3.11  shows  the  high  pressure  result  compared  with  atmos¬ 
pheric  pressure  data  obtained  by  other  investigators.  It  is  apparent  that 
if  the  pressure  does  affect  for  this  system  at  the  concentration  and 
temperature  considered,  the  effect  must  be  very  slight. 

This  represents  the  first  application  of  the  Agar-Turner  dual 
electrode  type  binary  Soret  cell  for  measurements  at  salt  concentrations 
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TABLE  3.1 


BINARY  THERMAL  DIFFUSION  EXPERIMENT 
AT  1000  BARS  PRESSURE 

EXPERIMENT  SDB-6 


Top  Electrode  Pair: 


t 

(t/0) 

exp- (t/0) 

R 

(«-R) 

ln|R-ig 

(minutes) 

(ohms) 

(ohms) 

9. 1 

.0502 

0.95103 

615.616 

3.056 

1.1171 

15.6 

.0861 

0.91750 

613.031 

0.471 

-0.7529 

25.6 

.  1413 

0.86822 

611.687 

-0.873 

-0.1358 

32.9 

.1816 

0.83393 

611.371 

-1.189 

0.1731 

38.7 

.2136 

0.80767 

611.254 

-1.306 

0.2670 

41.1 

.2268 

0.79708 

611.189 

-1.371 

0.3155 

47.9 

.2643 

0.  76774 

611.473 

-1.087 

0.0834 

61.5 

.3394 

0. 71219 

611.552 

-1.008 

0.0080 

69.2 

.3819 

0.68256 

611.545 

-1.015 

0.0149 

87.0 

.4801 

0.61872 

611.570 

-0.990 

-0.0100 

101.3 

.5591 

0.57172 

611.694 

-0.866 

-0. 1439 

115.3 

.6363 

0.52924 

611.770 

-0.790 

-0.2357 

130.2 

.7185 

0.48748 

611.871 

-0.689 

-0.3725 

140.5 

.7754 

0.46051 

611.819 

-0.741 

-0.2998 

155.6 

.8587 

0.42371 

612.000 

-0.560 

-0.5798 

175.3 

.9674 

0.38006 

612.038 

-0.522 

-0.6501 

192.6 

1.0629 

0.34545 

612.118 

-0.442 

-0.8164 

222.6 

1.2285 

0.29273 

612.110 

-0.450 

0.7985 

1121.4 

6. 1887 

0.00203 

615.628 

3.068 

1.  1210 

1127.4 

6.2219 

0.00200 

616. 140 

3.580 

1.2  754 

1138.3 

6.2820 

0.00190 

615.727 

3.167 

1.  1528 

Cond  it  ions : 

T 

27.253°C  :  AT  = 

9. 15°C 

YFeed 

3.57.  NaCl  :  B  = 

0.92  (see  reference  78,pp.3 

9 

181.20  minutes 

(x/a)  * 

0.81448  (print); 

0.85294  (base 

reversed) 

Derived  Results 

:  R 

(X 

612.56  ohms 

R 

o 

610.66  ohms 

°  12  * 

1.17  x  10"  3  °K'1 

;  1.048  x  10"3 

°K~  1 
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FIGURE  3.9.  EXTRAPOLATION  OF  THE  TOP  ELECTRODE  PAIR  RESISTANCE 

FOR  R  -  EXPERIMENT  SDB-6 
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FIGURE  3.10.  EXTRAPOLATION  OF  THE  TOP  ELECTRODE  PAIR  RESISTANCE 

FOR  R  -  EXPERIMENT  SDB-6  P11041 
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PRESSURE 


greater  than  0.05  molal  and  at  pressures  greater  than  atmospheric. 

Future  measurements  should  be  made  with  better  air  bath 
(surrounding  the  Soret  bomb)  control,  and  with  a  better  electrode  seal 
(the  bottom  electrode  pair  could  not  be  used  because  of  an  imperfect  seal 
between  the  epoxy  and  the  Monel). 

The  limitations  of  the  binary  Soret  cells  used  by  other  investi¬ 
gators  (90)  (125)  for  high  pressure  measurements  with  organic  liquids 
are  present  in  the  Soret  bomb  design  employed  for  this  study. 

3.4  EXTENSION  TO  MULTICOMPONENT  SYSTEMS 

Using  the  approach  of  Vinograd  and  McBain  (119),  Guthorie,  et  al. 
(122)  developed  a  relatively  simple  expression  for  the  ionic  Soret  coefficient 
in  a  dilute  ideal  multicomponent  ionic  system.  The  ionic  Soret  coefficient 
ia  defined  as. 


a 


(19) 


for  ion  i  in  the  system. 

A  series  of  atmospheric  pressure  experiments  were  performed 
using  different  multicomponent  Soret  cells  suggested  by  Hershey  and  Prados 
(123)  and  Berkau  and  Fisher  (124)  with  the  NaCl-I^O  system.  The  data  and 
auxiliary  information  are  reported  in  Table  3.2. 

From  these  experiments,  which  were  designed  to  determine  the 
Soret  cell  design  most  appropriate  for  multicomponent  aqueous  electrolyte 
solutions  such  as  sea  water,  it  was  found  that  metal  cell  parts  (such  as 
end  plates,  needles  sticking  in  from  the  side,  etc.)  should  not  contact 
the  multicomponent  salt  solution. 
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Since  for  binary  systems  (123)  (species  2  denotes  salt), 
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for  the  cell  design  w£.th  sample  removal  near  the  heated  and  cooled  plates, 
and  since  (124) 


CT 
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for  the  cell  design  with  a  center  dividing  membrane  (or  porous  barrier  to 
bulk  streaming),  it  is  possible  to  quantitatively  check  the  various  design 
alternatives  using  the  binary  Soret  coefficient  data  known  at  1  atm.  With 
Oj2  =  1.2  x  10  3  (°K  * )  (see  Figure  3.11), 


Y2  -  Y2  ^  +  0.04 
B  T 


and , 

Y2  -  Y-  ~  +  0.02 
B 

would  be  the  expected  concentration  differences  (in  weight  percent)  for  the 
end  sampling  and  center  dividing  type  multicomponent  Soret  cells,  respectively. 
Only  the  Teflon  wall,  center  dividing  type  cell  (without  metal  parts  exposed 
to  the  salt  solution)  yielded  measured  salt  separations  comparable  to  those 
expected  for  this  type  of  cell  (see  Experiments  MCSC-18  and  MCSC-19). 

Since  there  is  always  some  question  as  to  the  effect  of  a  cello¬ 
phane  diaphragm  on  ion  transport  rates,  a  porous  glass  diaphragm  would 
be  a  more  appropriate  choice. 

The  most  desirable  multicomponent  Soret  cell  for  high  pressure 
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(or  atmospheric  pressure)  measurements  is  shown  in  Figure  3.12.  It 
would  be  contained  in  a  cavity  similar  in  configuration  to  that  used  for 
the  binary  Soret  bomb  (see  Figure  3.6).  Oil  and  mercury  would  transmit 
the  pressure  to  the  solution  in  the  Teflon  cell  via  a  freely  moving 
o-ring  sealed  Teflon  piston. 

The  cold  chamber  would  have  a  much  greater  capacity  than  the 
hot  chamber  in  order  to  minimize  the  effect  of  solution  transport  through 
the  membrane  into  the  cold  chamber  during  depressurization.  No  difficulty 
would  be  encountered  during  pressurization  since  the  solution  in  both 
chambers  would  be  of  uniform  concentration. 
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FIGURE  3.12.  SORET  CELL  FOR  Ml’LTl COMPONENT  SALT  SOLUTIONS  AT  "1CH  PRESSURE 
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SECTION  4 


VISCOSITY  AND  THERMAL  CONDUCTIVITY 
OF  MULTICOMPONENT  SALT  SOLUTIONS 


4 . 1  INTRODUCTION 

In  connection  with  the  program  of  experimental  measurements  of 
the  ordinary  diffusion  coefficient  and  Soret  coefficient  in  binary  and 
multicomponent  salt  solutions,  a  review  and  analysis  of  existing  literature 
on  the  molecular  viscosity  and  thermal  conductivity  of  such  electrolytes 
recently  has  been  carried  out.  It  was  desired  to  relate,  where  possible, 
the  available  information  on  these  other  two  basic  transport  properties  to 
the  mass  diffusional  properties.  Both  experimental  viscosity  and  thermal 
conductivity  data,  and  methods  for  correlating  and  predicting  these  two 
properties  were  included  in  this  study. 

As  has  been  pointed  out  earlier  in  the  literature  (34)  (47),  the 
available  experimental  data  for  sea  water  itself  are  meager  with  regard 
to  viscosity,  and  absolutely  lacking  in  the  case  of  thermal  conductivity. 
Corresponding  data  for  other  multicomponent  salt  solutions  likewise  are 
sparse,  and  apparently  have  been  obtained  only  in  a  few  special  cases  of 
practical  interest  or  where  unusual  effects  were  anticipated.  Considerable 
experimental  data  on  the  viscosity  of  simple  binary  salt  solutions  are 
available,  however,  covering  ranges  of  concentration  and  temperature,  and, 
to  a  lesser  extent,  of  pressure.  Thermal  conductivity  data  for  binary 


electrolytes  are  somewhat  more  limited,  and  are  generally  reported  only 
for  various  concentrations. 

The  primary  objective  of  this  review  has  been  the  correlation 
and  prediction  of  the  viscosity  and  thermal  conductivity  of  multicomponent 
salt  solutions  as  functions  of  concentrations,  temperature,  and  pressure. 
Emphasis  here  has  been  given  to  the  study  of  sea  water  compositions  in 
the  temperature  range  0°C-30°C  and  the  pressure  range  1-1000  atm.  A 
secondary  and  continuing  objective  has  been, and  is, the  development  of 
theoretical  bases  for  relating  molecular  viscosity  and  thermal  conductivity 
to  the  mass  diffusional  transport  properties  being  studied  in  this  investi¬ 
gation. 

The  discussion  to  follow  will  include  a  preliminary  assessment 
of  sea  water  as  a  multicomponent  salt  solution,  a  review  and  analysis  of 
the  available  data  and  predictions  of  viscosity  and  thermal  conductivity 
for  those  systems  which  have  been  studied  in  detail,  and  an  outline  of 
theoretical  relationships  between  transport  properties  which  might  be 
useful  in  correlating  measurements  of  the  ordinary  diffusion  coefficient 
and  the  Soret  coefficient. 

4.2  COMPOSITION  OF  SEA  WATER 

Inasmuch  as  the  emphasis  on  multicomponent  salt  solutions  arises 
solely  from  the  eventual  goal  of  predicting  or  measuring  the  transport 
properties  of  sea  water,  a  brief  consideration  of  the  ionic  composition 
of  sea  water  should  be  appropriate  in  the  framework  of  the  discussion  of 
simpler  salt  systems  in  sections  to  follow.  Defant  (7)  has  given  typical 
compositions  of  sea  water  samples  having  salinities*  ranging  from  33  gm/kg 

» - 

*  "Salinity"  as  conventionally  defined  may  be  approximately  taken  to  be 
the  weight  of  dissolved  salts  (in  grams)  per  kilogram  of  sea  water. 
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up  to  38  gm/kg;  an  average  ionic  composition,  corrected  to  a  "standard" 
salinity  of  35  gm/kg,  is  shown  in  Table 4.1.  In  this  table,  the  ions  H+, 

HCO^,  and  H^BO^  are  listed  only  as  nominal  constituents,  inasmuch  as 
their  true  concentrations  (excluding  undissociated  CC^  and  H^BO^)  will 
depend  on  the  pH,  which  is  variable  from  point  to  point  in  the  oceans. 

Another  manner  of  considering  the  composition  of  sea  water 
might  be  in  terms  of  its  component  salts,  or,  rather,  in  terms  of  the 
various  binary  salt  solutions  which  would  have  to  be  mixed  together  in 
order  to  synthesize  sea  water.  One  possible  breakdown  along  these  lines 
is  shown  (for  a  salinity  of  35  gm/kg)  in  Table  4.2,  in  which  priority  has 
been  given  to  the  logical  "major"  binary  salt  components,  and  in  which 
the  total  ion  concentration  in  each  individual  solution  has  been  fixed 
at  the  sea  water  average  of  ~1 . 156M  for  this  salinity.  It  is  seen  that 
a  few  of  these  component  solutions,  such  as  "CaSO^"  and  are 

obviously  hypothetical  at  the  stated  concentrations  in  view  of  solubility 
limitations.  Nevertheless,  at  least  the  first  five  salts  (NaCl,  MgCl^. 
MgSO^,  CaSO^,  and  ^SO^)  are  predominantly  formed  upon  evaporation  of 
sea  water  to  dryness,  and  have  been  so  identified  in  the  solid  crystalline 
phase. 

It  is  apparent,  from  either  Table  4.1  or  4.2,  that  sodium 
chloride  is  by  far  the  predominating  binary  salt  component  of  sea  water. 

From  Table  4.2,  it  is  seen  that  the  sea  water  salt  system  is  approximately 
83  weight  percent  NaCl  and  11  weight  percent  MgC^,  so  that  these  two  salt 
components  total  nearly  94  percent  by  weight  of  the  entire  system.  Accord¬ 
ingly,  sea  water  itself  is  not  truly  a  very  good  example  of  a  complex 
multicomponent  salt  solution,  inasmuch  as  its  properties,  to  a  first 
approximation,  might  be  taken  equal  to  those  of  a  sodium  chloride  solution 
of  the  same  concentration,  and  furthermore,  to  a  second  and  somewhat  better 
approximation,  might  be  taken  equal  to  those  of  a  mixture  containing  a  major 
proportion  of  sodium  chloride  and  a  minor  proportion  of  magnesium  chloride. 
The  appropriate  validity  of  these  approximations  will  be  demonstrated  in 
the  discussions  to  follow. 


TABLE  4 . 1 


AVERAGE  IONIC  COMPOSITION  OF  SEA  WATER 
(S  -  Salinity  -  35  gm/kg) 


Molality  (~  Molarity, 
Ions  in  dilute  solution) 


Na+ 

0.4797 

K+ 

0.0103 

Mg** 

0.0553 

Ca^ 

0.0108 

Sr4* 

0.0002 

(H+) 

0.0035 

Cl" 

0.5641 

Br" 

0.0008 

so" 

4 

0.0291 

(HCO") 

0.0025 

(h2bo") 

0.0004 
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TABLE  A .  2 

HYPOTHETICAL  COMPOSITION  OF  SEA  WATER 
IN  TERMS  OF  COMPONENT  SALT  SOLUTIONS 

(S  «  Salinity  ■  35  gm/kg) 


Solution  Wt.  7.  of  Sea  Water 


0.578 

M 

NaCl 

82.72 

0.386 

M 

MgCl2 

10.97 

0.578 

M 

MgS04 

2.33 

0.578 

M 

"CaS04" 

1.94 

0.386 

M 

k2so4 

1.38 

0.578 

M 

"h2co3" 

0.43 

0.578 

M 

HBr 

0.14 

0.578 

M 

H3BO3 

0.04 

0.386 

M 

,,Sr(H2B03)2" 

0.03 

0.578 

M 

"SrS04" 

0.02 
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4.3  VISCOSITY  OF  ELECTROLYTE  SOLUTIONS 
4.3a  Viscosity  Data 

(1)  Sea  Water 

The  available,  published  data  on  the  viscosity  of  sea  water 
compositions  consist  of  the  original  observations  of  Krummel  (26)  and  the 
more  recent  values  of  Miyake  and  Koizumi  (33)  ,  Darmois  and  Darmois  (6) , 
and  Korf  and  Zakharova  (24).  The  first  two  of  these  four  cover  the  con¬ 
centration  range  from  zero  up  to  ~40  gm/kg  salinity  and  the  temperature 
range  0°C  to  30°C;  the  third  consists  only  of  a  series  of  observations  on 
undiluted  sea  water  at  28°C,  and  the  fourth  covers  a  higher  concentration 
range  of  sea  water  brines  and  the  temperature  range  -10°C  to  20°C.  In 
Figure 4.1  these  data  (with  the  exception  of  reference24,  which  could  not 
be  obtained  as  of  the  date  of  this  writing)  are  plotted  in  the  form  of 
"reduced  viscosity"  ur  versus  weight  percent  salts  in  solution  (Sa^^"^t:^)  , 
where  is  defined  as: 


in  which  iir  is  reduced  viscosity, 

is  measured  viscosity  of  solution  at  given  concentration 
and  temperature,  and 

M-w  is  measured  viscosity  of  pure  water  at  the  same 
temperature . 

Reduced  viscosity  is  a  more  sensitive  and,  at  the  same  time,  more  universal 
basis  for  comparing  data  of  several  investigators  than  would  be  the  measured 
absolute  viscosity  in  centipoises. 

The  lines  shown  in  Figure  4.1  represent  the  approximate  "best  fits" 
for  5°  temperature  increments  in  the  range  0°C-30°C.  Within  the  scatter 
of  the  data,  these  lines  of  reduced  viscosity  versus  concentration  were 
taken  to  be  linear  and  (with  the  exception  of  25°C  and  30°C  temperatures) 
were  drawn  so  as  to  pass  through  the  origin  (nr  *1).  In  subsequent  dis¬ 
cussion  it  will  be  shown  that  some  curvature  at  low  concentrations  would 
be  predicted  theoretically. 
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As  far  as  could  be  discerned  from  e  survey  of  the  literature, 
the  viscosity  of  sea  water  solutions  apparently  never  has  been  measured  as 
a  function  of  pressure,  or  for  temperatures  higher  than  30°C. 

(2)  Other  Multicomponent  Systems 

A  limited  number  of  investigations  of  the  viscosities  of 

other  multicomponent  salt  solutions  have  been  reported  in  the  literature. 

In  a  recent  study  of  viscosities  of  the  systems  CdClj-MCl-l^O,  where 

M  «  H+,  Li+,  Na+,  K+,  and  ,  Galinker,  Tyagai,  and  Fenerli  (15)  found 

that  in  every  case  the  viscosities  of  mixtures  were  lower  than  might  have 

(  |  | 

been  predicted  from  simple  additivity,  owing  to  the  formation  of  MCd 
complex  ions.  In  an  earlier  investigation,  Tollert  (45)  had  observed  a 
similar  interaction  in  the  case  of  MgCl_-KCl-H_0,  indicating  the  presence 

j  |  |  ^ 

of  the  complex  ion  KMg  .  Based  upon  these  findings,  some  interaction 

influence  in  the  case  of  sea  water  might  be  expected  by  virtue  of  formation 

|  |  ) 

of  the  complex  cation  NaMg  ,  but  any  such  effect  would  be  relatively  minor 

in  view  of  the  nearly  10:1  non-stoichiometric ,  disproportionate  excess  of 

„  +  „  -H- 

Na  over  Mg 

(3)  Binary  Salt  Solutions 

As  might  be  expected,  the  great  bulk  of  experimental  data 
on  the  viscosity  of  aqueous  electrolyte  solutions  pertain  to  simple,  binary 
salt  solutions.  Of  primary  interest  here  will  be  such  data  on  sodium 
chloride  solutions  at  various  concentrations  and  temperatures.  In  Figure  4.2 
is  shown  a  comparison  between  the  reduced  viscosities  of  sea  water  solutions 
and  sodium  chloride  solutions  at  concentrations  from  zero  to  20  weight  per¬ 
cent,  and  for  temperatures  in  the  range  0°C-60°C,  where  the  sodium  chloride 
lines  represent  the  smoothed  and  averaged  data  of  Sergeevich ,Zhuze ,  and 
Chestnov  (41,  48)  and  of  Suryanarayana  and  Venkatesan  (42).  The  latter 
data  actually  extend  up  to  saturation  (  ~26-27  weight  percent  NaCl) ,  and 
will  be  shown  for  this  range  in  a  subsequent  plot;  a  direct  comparison  with 
sea  water  in  the  vicinity  of  "saturation"  would  not  be  meaningful,  however, 
inasmuch  as  sea  water,  as  a  multicomponent  mixture,  does  not  have  a 
unambiguous ly-def inable  "solubility  limit". 
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AX  ISO  DS  ia  asonasa  ,Jrf 


FIGURE  4.2.  COMPARISON  OF  SMOOTHED  EXPERIMENTAL  VISCOSITY  DATA  FOR  SEA 
WATER  AND  AQUEOUS  SODIUM  CHLORIDE  SOLUTIONS  AT  ATMOSPHERIC 
PRESSURE 


It  is  seen  from  Figure  4.2  that  sea  water  solutions  at  any  given 
weight  percent  concentration  and  temperature  are  more  viscous  than  the 
corresponding  sodium  chloride  solutions.  In  the  concentration  region  of 
normal  sea  water  (  ~35  mg/kg  salinity  or  3.5  weight  percent),  the  viscosity 
ratio  of  sea  water  to  sodium  chloride  is  approximately  1.020  at  all  tem¬ 
peratures.  As  will  be  shown  later,  this  difference  is  that  which  would 
be  expected  from  the  contributions  of  the  minor  salt  ions  in  sea  water. 

(4)  Pressure  Effect 

The  viscosity  data  presented  up  to  this  point  have  all 
represented  values  at  essentially  atmospheric  pressure.  In  considering 
the  effect  of  pressure  on  viscosity,  it  is  instructive  to  first  consider 
pure  water,  for  which  reasonably  extensive  data  are  available.  In  Figure  4.3 

Up 

is  shown  a  plot  of  "relative  viscosity"  (— rr^ )  versus  pressure  for  pure 

h-q 

water  at  temperatures  in  the  range  0°C  to  75°C,  where  M-p  is  the  viscosity 
at  pressure  "p"  and  y.Q  is  the  viscosity  at  zero  pressure  (essentially  that 
at  1  atm).  The  curves  shown  represent  averages,  at  least  in  the  lower  - 
pressure  range,  between  the  values  of  Bridgman  (4)  and  the  earlier  values 
of  Cohen  (5) .  It  is  seen  that  viscosity  rises  with  pressure  at  the  higher 
temperatures,  whereas  at  the  lower  temperature  (<~  30°C)  there  is  a  small 
initial  decrease  in  viscosity,  followed  by  a  subsequent  rise. 

For  oceanographic  purposes  only  the  pressure  range  1-1000  atm 
is  of  general  interest.  Figure 4.4  expands  this  relatively  small  region  of 
Figure  4, 3, for  which  the  relative  viscosity  of  pure  water  varies  only  from 
0.92  to  1.08  over  the  entire  temperature  range  0°C  to  75°C.  In  Figure  4.4  and  4.5 
are  shown,  for  comparison  with  the  water  data,  the  relative  viscosities  of 
sodium  chloride  solutions  in  the  concentration  range  0-20  weight  percent 
and  the  temperature  range  5°C  to  40°C,  representing  the  smoothed  data  of 
Sergeevich,  Zhuze,  and  Chestnov  (41,  48)  for  pressures  up  to  500  atm. 

(It  is  understood,  in  the  definition  of  "relative  viscosity"  for  salt 
solutions,  that  "pip"  and  "n0"  both  refer  to  solutions  of  a  given  concen¬ 
tration,  rather  than  to  pure  water.)  It  is  seen  that  the  effect  of  adding 
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(AVERAGED  DATA  OF  BRIDGMAN  AND  OF  COHEN) 


VARIOUS  TEMPERATURES 


(SMOOTHED  DATA  OF  ZHUZE,  SERGEEVICH,  AND  CHESTNOV) 


FIGURE  4.4.  RELATIVE  VISCOSITIES  OF  AQUEOUS  SODIUM  CHLORIDE  SOLUTIONS  IN  THE  PRESSURE 
RANGE  0  -  1000  ATM,  TEMPERATURE  RANGE  5°C  -  40°C 


(SMOOTHED  DATA  OF  ZHUZE,  SERGEEVICH,  AND  CHESTNOV) 


FIGURE  4.5.  RELATIVE  VISCOSITIES  OF  AQUEOUS  SODIUM  CHLORIDE  AND 
PURE  WATER  IN  THE  PRESSURE  RANGE  0-1000  ATM, 
TEMPERATURE  RANGE  5°C  -  40°C 


salt  to  water  is  to  make  the  pressure  effect  more  positive  in  every  in¬ 
stance.  For  example,  at  a  temperature  of  20°C,  the  viscosity  of  pure  water 
at  500  atm  pressure  is  1.4  percent  less  than  at  1  atm,  whereas  that  of 
3.5  weight  percent  sodium  chloride  solution  at  500  atm  pressure  is  1.4  per¬ 
cent  greater  than  at  1  atm. 

b)  Viscosity  Correlations 

Numerous  methods  for  correlating  the  viscosities  of  electrolyte 
solutions  versus  concentration,  temperature  and  pressure  have  been  proposed, 
some  having  at  least  partial  theoretical  justification,  and  others  being 
no  more  than  strictly  empirical.  Of  these  various  correlations,  the  more- 
accepted  ones  recently  have  been  reviewed  in  a  survey  article  by  Kozlowska  (25) . 
The  method  consistently  receiving  the  greatest  attention,  at  least  in  the  cor¬ 
relation  of  viscosity  versus  concentration,  has  been  that  based  upon  the 
original  equation  of  Jones  and  Dole  (18) ,  which,  in  its  more  generalized 
form,  is  now  written  as  follows: 

H  -  1  +  BC  +  DC2  (4-2) 

r 

in  which  M-r  is  reduced  viscosity, 

C  is  concentration,  moles/liter, 

A,  B  and  D  are  coefficients,  having  the  units  (moles/liter) 

-1  -2 

(moles/liter)  ,  and  (moles/liter)  ,  respectively. 

The  "A"  coefficient  in  this  equation  can  be  calculated  theoretically  for 
simple  systems,  while  any  determination  of  the  "B"  and  "D"  coefficients 
is  at  least  partly  empirical,  and  must  be  based  on  experimental  values. 

While  other  correlations  of  viscosity  versus  concentration  have  been  sug¬ 
gested,  they  generally  have  been  found  to  be  less  satisfactory. 

The  variation  of  viscosity  of  electrolyte  solutions  with  temper¬ 
ature  normally  is  correlated  by  means  of  an  Arrhenius -type  relationship, 
in  which  viscosity,  in  effect,  is  considered  as  a  reciprocal  ’’rate  constant" 
for  viscous  flow.  The  Arrhenius  activation  energy  thus  obtained  for  the 
solution,  however,  generally  has  only  been  related  semiempirically  to  the 
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electrolyte  concentration,  as,  for  example,  in  terms  of  the  percentage  of 
saturation  (27,  44).  The  relatively  meager  data  available  on  the  vis¬ 
cosities  of  electrolyte  solutions  at  elevated  pressures  has  been  corre¬ 
lated  with  the  changes  of  conductivity  resulting  from  increases  in 
pressure  (35) . 

Most  correlations  or  predictions  of  viscosity  have  been  applied 
to  the  relatively  less-complicated  case  of  simple  binary  electrolytes. 

Those  few  correlations  valid  for  mixtures  generally  have  been  based  upon 
some  sort  of  additivity  principle,  in  which  experimental  knowledge  of  the 
contributions  of  each  individual  ion  or  salt  component  is  required  as  input. 

The  various  methods  for  correlating  orjpredicting  the  viscosities 
of  aqueous  electrolyte  solutions  will  be  discussed  in  further  detail,  and 
inferences  will  be  drawn  where  appropriate  regarding  their  applicability 
toward  the  sea -water  system. 

(1)  The  Jones -Dole  Correlation 

The  "complete"  Jones -Dole  correlation  expressed  by  equation 
4.2  can  be  considered  in  somewhat  simpler  form  in  regions  of  either  very 
low  or  ”ery  high  concentration.  The  limiting  law  first  proposed  by  Falken- 
hagen  (10)  for  the  case  of  extreme  dilution  is  the  following: 

M>  ■  1  +  (4.3) 

r 

From  considerations  of  Debye-Huckel  theory,  in  which  the  ions  are  considered 
to  act  as  point  charges,  the  coefficient  "A"  can  be  predicted  theoretically, 
either  for  binary  electrolytes,  or,  by  suitable  definition  of  the  "ionic 
strength"  for  mixtures,  The  positive  coefficient  "A"  generally  has  a 
numerical  value  in  the  vicinity  of  0.006  (moles/liter)  ^  for  monovalent 
electrolytes,  and  in  the  range  0.01-0.03  (moles/liter) for  electrolytes 
containing  polyvalent  ions  (25).  Thus,  for  the  concentration  region 
(<0.01M)  where  equation  4.3  by  itself  truly  does  represent  the  sole  effect 
of  concentration  on  viscosity,  the  maximum  measurable  increase  in  viscosity 
over  that  of  pure  water  is  only  a  small  fraction  of  one  percent. 
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For  electrolyte  concentrations  in  the  range  0.01M  <  C  <  1.0M 
the  "classical"  form  of  the  Jones-Dole  equation  is  a  reasonably  good 
approximation: 

M,  =  1  +  A\£~+  BC  (4.4) 


At  concentrations  much  higher  than  this,  however,  a  small  quadratic  term 
is  required,  as  indicated  by  equation4.2  for  the  general  case.  One  mis¬ 
conception  which  is  sometimes  taken  for  granted  is  that  both  of  the  non¬ 
linear  terms  AyC” and  DC^  can  be  completely  neglected  throughout  an  inter¬ 
mediate  concentration  region  (i.e.,  from  0.1M  to  1.0M),  so  that  the 
viscosity  equation  could  be  written  simply  as: 


1  +  BC 


(4.5) 


To  be  sure,  viscosity-concentration  relationships  may  be  essentially 
linear  in  this  region,  but  the  lines  thus  obtained  do  not  extrapolate 
back  to  the  origin  (p^  ■  1).  More  correctly,  the  quasi-linear  relation 
at  intermediate  concentration  should  be  written  as: 


*  (1  +  X)  +  (B  +  Y)C 


(4.6) 


in  which  "X"  and  "Y"  represent  small,  but  not  completely  negligible  con¬ 
tributions  from  the  non-linear  terms.  Indeed,  if  equation4.6  is  made  to 
satisfy  exactly  both  value  and  derivative  at  the  inflection  point  of  the 
true  equation  4.2: 

i  /.  /  i  i/-i 

(4.7) 


X  -  (~)(A)4/3(D)'1/3 


(f)  (A)2/3(D) 1/3 


(4.8) 


Unlike  the  coefficients  "A"  and  "D",  both  of  which  are  always 
positive,  the  coefficient  "B"  in  the  Jones-Dole  equation  can,  in  relatively 
rare  instances,  take  on  a  negative  value.  In  general,  however,  the  co¬ 
efficient  "B"  is  positive,  and  of  sufficient  order  of  magnitude  to  account  for  most 
of  the  observed  Increase  in  the  viscosity  Of  electrolytes  over  that  of  pure 
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water.  Ion  hydration  (particularly  of  cations)  has  been  assumed  to 
account  for  the  major  positive  contfcibution  to  the  "B"  coefficient, 
whereas  breaking  of  the  structure  of  the  adjacent  water  (particularly 
by  anions)  may  lead  to  what  is  observed  as  a  “negative  viscosity" 
effect  (19,  20).  Both  of  these  influences  may  be  classed  as  "ion-solvent" 
interaction,  in  contrast  to  the  "ion-ion"  interaction  which  determines 
the  value  of  the  "A"  coefficient. 

Several  correlations  of  the  Jones  and  Dole  "B"  coefficient 
versus  other  measurable  or  calculable  properties  of  electrolyte  solu¬ 
tions  have  been  proposed.  Bramhall  (3),  for  example,  has  proposed  an 
empirical  correlation  between  viscosity  and  density,  in  which  relative 
increase  (or  decrease)  in  viscosity  coefficient  "B"  is  found  to  be  pro¬ 
portional  to  the  percentage  contraction  (or  expansion)  of  the  water  volume 
fraction.  Fisher  and  Zaitseva  (13),  on  the  other  hand,  lave  derived  a 
theoretical  procedure,  based  upon  a  hydrodynamic  theory  of  ion  motion, 
for  estimating  the  contribution  of  ion  hydration  alone  uo  the  overall 
"B"  coefficient;  their  method,  however,  does  not  account  for  the  Coulomb 
ion-ion  interaction  contribution  in  the  concentration  region  where  the 
Debye -Hu’ckel  formalization  begins  to  break  down. 

Perhaps  the  most  widely  used,  albeit  strictly  empirical , method 

for  estimating  values  of  the  "B"  coefficient  for  binary  electrolytes  has 

been  that  based  upon  the  additivity  of  individual  ,,B^on"  contributions. 

It  has  been  found,  for  example,  that  the  viscosity  of  aqueous  KjSO^  can 

be  estimated  reasonably  well  from  the  additivity  principle  and  the  known 

viscosity  behavior  of  aqueous  KC1,  NaCl,  and  Na^SO^.  Inasmuch  as  no  ion 

in  solution  can  exist  by  itself,  the  assignment  of  individual  "B,  " 

ion 

values  must  be  somewhat  arbitrary.  Kaminsky  (19,  20)  has  proposed  such 
an  assignment  based  upon  the  starting  assumption  that  B^  =  E_^-  = 
at  all  temperatures. 

A  natural  temptation  with  the  "B.  "  method  would  be  to  assume 

ion 

that  the  additivity  principle  (based  on  ionic  molarities)  holds  true  for 
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multicomponent  as  well  as  for  binary  electrolytes.  As  will  be  pointed  out 
later,  this  assumption  is  probably  an  oversimplification;  but  nonetheless, 
for  a  mixture  such  as  sea  water  (which  is  dominated  by  the  two  major  ions 
Na+  and  Cl"),  it  may  yield  surprisingly  accurate  predictions.  A  compari¬ 
son  between  reduced  viscosity  values  for  sea  water  solutions  (a)  pre¬ 
dicted  from  this  assumption  and  (b)  representing  averaged  experimental 
data  accordingly  is  shown  in  Figure  4.6  for  the  concentration  range 
0-4  percent,  and  the  temperature  range  0  -  35°C.  By  comparison  of  Figure 
4.6  with  Figure  4.1,  it  is  seen  that  the  predictions  generally  agree 
with  the  experimental  data  within  the  limits  of  uncertainty  of  the  data 
itself . 

(2)  Other  Viscosity-Concentration  Correlations 

Several  authors  (25,  35)  concerned  with  electrolyte  vis¬ 
cosity  have  chosen  to  apply  the  classical  Jones-Dole  equation  in  its 
"inverse"  form,  i.e.,  in  terms  of  reduced  fluidity  rather  than  reduced 
viscosity: 

—  =  — ~  ~  1  -  AVC"  -  B*C  -  D'C2  (4.9) 

in  which  B'  and  D'  are  new  coefficients  having  the  units  (moles/liter)  * 

-2 

and  (moles/liter)  ,  respectively.  In  the  region  of  extreme  dilution  where 

Falkenhagen  theory  is  applicable,  the  relative  changes  in  viscosity  are  so 

small  that  the  coefficient  "A"  has  essentially  the  same  value  whether 

equation  4. 2  or  equation  4.9  is  used.  The  empirically  determined  "B"'  and 

"D'"  coefficients,  on  the  other  hand,  are  somewhat  different  numerically 

3/2 

than  their  "B"  and  "D"  counterparts.  Indeed,  if  the  missing  "C  "  term 
and  corresponding  higher  terms  are  added  to  create  a  homologous  series, 
the  new  coefficients  of  the  "inverse"  Jones-Dole  equation  then  can  be 
exactly  related  mathematically  to  the  coefficients  of  the  "conventional" 
Jones-Dole  equation. 

Another  type  of  viscosity-concentration  correlation  has  been  the 
empirical  "exponent"  method  proposed  by  Andrussow  (1,2).  In  this  correlation, 
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viscosity  and  concentration  are  related  exponentially  as  follows: 


(4.10) 


in  which  42  and  4^  ate  viscosities  at  concentrations  C2  and  C^,  respectively, 
and  where  the  exponent  "Y"  itself  is  correlated  by  means  of  a  "Jones -Dole" 
type  expansion: 


Y~  Yc  +  a(C2-C1)1/2  +  b^-C^)  +  c(C2-C1)3/2  +  .  (4.11) 

(units  of  a,  b,  c,  etc.  are  (moles/liter)  , 

(moles/liter)  *,  (moles/liter) ~^/2 f  etc.) 

Normally,  the  "reference"  concentration  "C^"  is  fixed  at  some  arbitrarily 
low  level  (M3.001M)  in  the  region  where  viscosity  can  be  described  by  the 
ideal  Falkenhagen  theory. 

The  advantage  claimed  by  Andrussow  for  the  "exponent"  method  was 
that  the  empirical  coefficients  "a",  "b",  "c",  etc.  were  much  less 

dependent  on  temperature  than  are  the  coefficients  "A",  "B",  and  "D"  of 
the  Jones-Dole  equation.  Offsetting  this  advantage,  however,  as  later 
pointed  out  by  Kaminsky  (21)  is  the  disadvantage  that  the  coefficients 
"Y0",  "a",  "b" ,  "c" ,  etc.  are  truly  constant  only  over  a  relatively  limited 
range  of  concentrations,  so  that  several  sets  of  coefficients  must  be 
specified  to  cover  the  entire  concentration  range  from  zero  up  to  saturation. 
Indeed,  if  equations  4.3,  4.10  and  4.11  are  combined  and  simplified  for  the  case 
of  very  small  values  of  "C^",  the  result  is: 


U2  -C 

—  -  1  +  1 


{c^[Yo  +  *^2  +  bC2  +  cC2 


3/2 


]«. 


12) 


r2  -w 

from  which  it  is  seen  that  any  direct  correspondence  with  the  general  form 
of  the  Jones-Dole  correlation  is  impossible,  owing  to  the  presence  of  the 


highly  non-linear 


factor. 
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One  final  method  of  correlating  the  viscosities  of  both  electro¬ 
lytes  and  non-electrolytes  as  functions  of  concentration,  especially  for 
concentrations  in  the  vicinity  of  saturation,  is  that  originally  proposed 
by  Suryanarayana  and  Venkatesan  (42,  43),  as  follows: 

(4.13) 

in  which  ps  is  viscosity  of  the  saturated  solution  at  a  given  temperature, 
or  and  3  are  constants  for  any  given  temperature,  and 

C  is  the  saturation  solubility  at  the  given  temperature, 
s 

Kume  and  Tanaka  (27,  44)  later  attempted  to  generalize  and  improve  this 
type  of  correlation  to  include  temperature  effects  as  well,  and  did  find 
that  for  most  aqueous  electrolytes  a  plot  of  "lnlj—  j"  versus  "  J  "  yielded 
points  on  a  single  smooth  (but  not-quite-linear)  cutve,  irrespective  of 
either  temperature  or  concentration,  rather  than  on  a  series  of  straight 
lines  for  different  temperatures,  as  would  be  predicted  from  equation  4.13 
They  were  unable,  however,  to  characterize  these  curves,  except  empirically 
for  each  electrolyte. 

(3)  Viscosity* Temperature  Correlations 

Most  correlations  of  electrolyte  viscosity  versus  temperature 
have  been  based  upon  equations  of  the  Arrhenius  type,  similar  to  those  used 
to  correlate  the  viscosities  of  pure  liquids: 

EvU./RI 

u  v  Ke  (4.14) 

in  which  pi  is  viscosity,  cp. 

K  is  pre-exponential  factor,  cp. 

Evis.  is  energy  of  activation  for  viscous  flow,  kcal/mole. 

T  is  temperature,  °K. 

R  -  0.001987  Kcal/mole-°K. 
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The  activation  energy  for  viscous  flow  MEv^a  ”  can  be  predicted  theoreti¬ 
cally.  in  the  case  of  non-associated  liquids,  by  means  of  the  theory  of 
rate  processes  originally  developed  by  Eyring  (8) . 

In  Figure 4.7,  Arrhenius -type  viscosity-temperature  correlations 
are  compared  for  pure  water,  normal  sea  water  (S  ■  35  gm/kg) ,  and  aqueous 
sodium  chloride  solutions  from  3.5  weight  percent  concentration  up  to 
saturation.  It  is  seen  that  the  electrolyte  viscosities  as  functions  of 
temperature  generally  parallel  that  of  pure  water,  so  that  any  effect  of 
concentration  on  the  total  activation  energy  for  viscous  flow  must  be 
relatively  minor.  It  is  also  apparent  that  the  activation  energy  "Ev^g  " 
for  water  alone  varies  rather  substantially  with  temperature,  ranging 
from~5.0  Kcal  at  0°C  down  to  ~3.0  Kcal  at  100°C.  The  variation,  however, 
is  continuous  throughout  the  entire  temperature  regime;  there  is  no  un¬ 
usual  effect  at  ^r  below  0°C  (as  is  the  case  with  electrolytic  conductivity). 

It  is  evident  from  equation  4.14  that  if  both  the  viscosity  of 
water  and  the  viscosity  of  aqueous  electrolytes  approximately  follow  the 
Arrhenius  law,  reduced  viscosity  "nr"  also  should  follow  this  law.  From 
Figures  4.1  and  4.2  it  can  be  seen  that  the  reduced  viscosities  of  sea  water 
and  of  aqueous  sodium  chloride  generally  increase  slightly  with  rising 
temperature,  so  that  the  total  activation  energies  "Ev^g  "  for  the  electro¬ 
lyte  solutions  are  slightly  less  than  that  of  pure  water.  However,  attempted 
Arrhenius  correlations  of  reduced  viscosities,  versus  temperature  showed  that 

"AE  "  (the  difference  in  activation  energy  between  the  electrolyte  and 
VIS* 

pure  water)  varied  too  widely  over  even  small  ranges  of  temperature  for 
this  type  of  correlation  to  be  of  much  value.  Kaminsky  (19,  20)  en¬ 
countered  similar  behavior  in  attempting  to  correlate  the  Jones-Dole 
individual  "B^  "  contributions  versus  temperature.  From  consideration  of 

equation  4.2  for  the  generalized  Jones-Dole  correlation,  it  is  of  course 
apparent  that  reduced  viscosity  "ur"  and  Jones-Dole  coefficient  "B"  cannot 
both  follow  the  Arrhenius  relationship  simultaneously. 
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FIGURE  4.7.  VISCOSITY-TEMPERATURE  VARIATION  FOR  PURE  WATER, 
SEA  WATER,  AND  AQUEOUS  SODIUM  CHLORIDE  SOLUTIONS 
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The  "degree  of  saturation"  correlations  of  viscosity  versus 
concentration  by  Kume  and  Tanaka  (27,  44)  have  already  been  mentioned. 

In  correlating  certain  electrolyte  viscosities  versus  temperature,  they 
found  that  a  better  Arrhenius -type  fit  could  be  obtained  if  viscosity  at 
a  constant  percentage  of  saturation  were  correlated,  rather  than  viscosity 
at  a  fixed  nominal  concentration  (weight  percent  or  molality).  The  re¬ 
sulting  equation  was: 

[E  (S  )] 
r  x 

- — 

(Hr)c  -  [Kr(Sx)]e  (4.15) 

cl " s* 

in  which  is  "degree  of  saturation",  and  where  Kf(Sx)  and  Er(Sx)  are 
empirical  functions  of  which  are  characteristic  for  any  electrolyte. 

Rao  and  Panicker  (37)  determined  the  viscosities  of  certain 
supersaturated  electrolyte  solutions,  and  found  that,  at  a  fixed  concen¬ 
tration,  viscosity  followed  a  decreasing  linear  relationship  with  rising 
temperature  for  some  electrolytes,  and  the  more  usual  exponential  or 
"Arrhenius"  relationship  with  temperature  for  others.  They  apparently 
did  not  observe  any  anomolous  behavior  of  viscosity  in  the  vicinity  of 
the  saturation  point  itself.  Fujita  (14),  on  the  other  hand,  found 
apparent  small  discontinuities  in  the  activation  energies  for  viscous  flow 
of  dilute  Na2C0^  solutions  at  the  two  temperatures  where  corresponding  dis¬ 
continuities  in  saturation  solubility  occur,  and  cited  similar  observed 
abnormalities  of  electrolytic  conductivity  for  dilute  Na2S0^  solutions.  It 
is  not  known  whether  any  such  anomolous  behavior  exists  in  the  case  of  a 
complex  mixture  such  as  sea  water. 

(4)  Viscosity-Pressure  Correlations 

Aside  from  theoretical  relationships  between  viscosity  and 
pressure  proposed  for  pure,  nonassociated  liquids  (4),  little  attention 
has  been  given  to  correlating  viscosity  with  pressure  for  aqueous  solutions 
of  electrolytes.  Podolsky  (35)  has  outlined  a  method  based  upon  experi¬ 
mental  knowledge  of  the  effect  of  pressure  on  electrical  conductance,  but 
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was  able  to  test  it  only  in  the  two  cases  of  aqueous  KC1  and  aqueous  NaCl. 

In  view  of  the  relatively  poor  agreement  between  data  and  theory  apparently 
obtained  in  the  latter  case,  this  correlation  may  not  be  of  too  much  value 
in  predicting  the  very  small  second-order  effect  of  electrolyte  concentration 

on  relative  viscosity  at  elevated  pressure. 

(5)  Viscosity  Correlations  for  Multicomponent  Electrolytes 

As  pointed  out  in  the  previously-cited  review  article  by 
Kozlowska  (25),  the  viscosities  of  multicomponent  electrolyte  solutions 
have  been  correlated  empirically  both  on  a  basis  of  ionic  strength  and 
on  a  basis  of  total  molarity  (or  molality).  Galinker,  Tyagai  and  Fenerli  (15) 
recently  proposed  some  alternative  "additivity"  relations  (for  mixtures  of 
non-polar  liquids  and  the  like)  which  might  be  valid  for  aqueous  electrolytes 
as  well: 


Additivity  of  viscosities: 


Vl  +  “2X2  + 


Additivity  of  cube  roots  of  viscosities: 
1/3  ~  1/3  1/3 

li  =  ^1  xi  +  ^2  *2  *  . 


Additivity  of  logarithms  of  viscosities: 
log  u  ~  x^og  u  +  x2  log  42  + 


Additivity  of  kinematic  fluidities: 


V 


+  x 


2  k 


(A. 16) 


(4.17) 

(4.18) 


(4.19) 


in  which  u  is  viscosity. 

x  is  mole  fraction  ,  and 
p  is  density  ; 

Subscripts  "1",  "2",  etc.  refer  to  component  1,  component  2,  etc. 

In  applying  these  four  relations  to  multicomponent  electrolyte 
solutions,  of  course,  the  "components"  would  be  understood  to  be  the 
minimum  necessary  number  of  binary  salt  solutions  (of  roughly  equal  con¬ 
centrations)  to  make  up  the  mixture,  rather  than  merely  the  individual 
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solid  crystalline  salts  plus  water.  Accuracy  moreover  could  be  improved 
by  replacing  absolute  viscosities  with  reduced  viscosity  differences, 


as  follows: 

Additivity  of  viscosities: 

(u  -1)  -  (u  -l)x.  +  (u  -l)x_  +  .  (A. 20) 

r  r  ^  i  z 

Additivity  of  cube  roots  of  viscosities: 

(nr1/3-l)  ■  (14. 1/3-l)x  +  (u  1/3-l)x  +  .  (4.21) 

r  1  1  r2  * 

Additivity  of  logarithms  of  viscosities: 

lo®10**r  *  xllo*10tlr1  +  x2log10,1r2  +  .  (4.22) 


Additivity  of  kinematic  fluidities: 


/M 

t 

U 

dL 

+  x2 

2  ' 


2' 


♦ 


(4.23) 


The  various  "additivity"  rules,  equations  4.20  through  4.23,  were  tested 
for  the  sea  water  system,  and  were  found  to  yield  essentially  equally  good 
predictions  of  viscosity.  Equation  4.20,  of  course,  will  be  recognized  to 
be  equivalent  to  the  additivity-of ""BjonM  method  used  in  calculating  the 
predictions  shown  in  Figure  4.6.  Galinker  et  al  .(15) ,  on  the  other  hand,  postu¬ 
lated  that  the  latter  three  relations  are  better  for  predicting  mixture 
viscosity  in  the  general  case.  In  view  of  the  relative  domination  of  the 
sea  water  system  by  the  two  ions  Na+  and  Cl  ,  linear  additivity  (equation  4.20) 
is  probably  as  good  an  assumption  as  would  ever  be  required  for  any  practical 
purpose . 

(6)  Sample  Predictions  of  Sea  Water  Viscosity 

Inasmuch  as  the  theoretical  or  empirical  methods  of  prediction 
just  discussed  appear  to  offer  little  or  no  advantage  over  what  experimental 
data  is  available,  the  recommended  method  embodies  use  only  of  experimental 
correlations . 
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(a)  Predict  the  viscosity  of  ordinary  sea  water  (S  »  35  gm/kg) 
at  5°C  and  500  atm. 

M-  of  ^0  at  5°C,  1  atm  =  1.519  cp  [Figure  4. 7] 

of  3.57.  sea  water  at  5°C,  1  atm  «  1.059  [Figure  4.1] 
U  n 

of  3.57.  NaCl  at  5  C  and  500  atm  -  0.976  [Figure  4.4] 
'o 

Mi  of  3.57.  sea  water  at  5°C  and  500  atm 
-  (1.519) (1.059) (0.976)  -  1.570  cp. 


(b) 


Predict  the  viscosity  of  107.  sea  water  brine 
li  of  ^0  at  60°C,  1  atm  ■  0.470  cp. 

(ur)8ea  ~ 

-z — r at  any  given  weight  concentration  ^ 

^r'NaCl 

of  107.  NaCl  at  60°C  and  1  atm  -  1.212 

4  of  107.  sea  brine  at  60 °C  and  1  atu 
-  (0.470) (1.020) (1.212)  -  0.581  cp. 


at  60°C  and  1  atm. 
[Figure  4.7] 

1.020 

[Figure  4.2] 
[Figure  4.2] 


4.4  THERMAL  CONDUCTIVITY  OF  ELECTROLYTE  SOLUTIONS 
a.  Thermal  Conductivity  Data 
(1)  Sea  Water 

No  published  data  apparently  are  available  on  the  thermal 
conductivity  of  sea  water.  Krummel  (26)  once  estimated  thermal  conductivity 
values  at  17.5°C  for  sea  water  salinities  from  zero  up  to  S  *  35  gm/kg  by 
making  the  rather  questionable  assumption  that  the  thermal  diffusivity  of 
sea  water  was  everywhere  equal  to  that  of  pure  water.  His  resulting  values 
are  plotted  in  Figure  4.8  for  comparison  with  presently  predicted  sea  water 
values  and  available  experimental  data  for  aqueous  NaCl. 
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(TEMPERATURE:  20°C) 


(2)  Binary  Salt  Solutions 

A  fair  amount  of  thermal  conductivity  data  for  binary 
electrolyte  solutions  Is  available,  although  most  of  It  pertains  only 
to  room-temperature  conditions.  In  general,  with  only  a  few  exceptions, 
thermal  conductivity  drops  with  increasing  electrolyte  concentration, 
although  the  relative  effect  is  much  smaller  than  is  the  case  for  vis¬ 
cosity.  Riedel  (39)  has  tabulated  thermal  conductivity  data  for  NaCl 
solutions  (along  with  many  other  salts,  acids,  and  bases)  at  concen¬ 
trations  from  zero  up  to  saturation  and  for  a  temperature,  in  general, 
of  20°C.  For  all  practical  purposes  the  change  (decrease)  in  thermal 
conductivity  for  most  salts  appears  to  be  linear  with  concentration 
(molarity)  throughout  the  entire  range. 

In  Figure  4. 8  the'^educed”  thermal  conductivity  function  (kf) 
for  sodium  chloride  is  plotted  versus  concentration  in  the  range  zero 
to  four  percent  by  weight,  and  with  "krM  defined,  analogously  to  "ur", 
as  follows: 

kr  -  rp  (4.24) 

w 

in  which  kf  is  reduced  thermal  conductivity 

k  is  thermal  conductivity  of  electrolyte,  cal/cm-sec-°C . 
k^  is  thermal  conductivity  of  pure  water,  cal/cm-sec-°C. 

It  is  seen  that,  at  the  3.3  weight  percent  concentration  level,  the  change 
in  thermal  conductivity  of  aqueous  NaCl  relative  to  that  of  pure  water 
amounts  to  only  0.6  percent.  In  subsequent  discussion,  a  comparison  will 
be  made  between  experimental  thermal  conductivity  values  for  NaCl  and 
predicted  thermal  conductivities  for  sea  water,  as  indicated  in  Figure  4.8. 

On  the  basis  of  thermal  conductivity  measurements  for  relatively 
concentrated  solutions  of  NaOH  and  KOH  in  the  temperature  range  0°C  to 
80°C,  Riedel  (38)  concluded  that  any  effect  of  temperature  on  reduced 
thermal  conductivity  "kr"  for  these  two  systems  was  negligible,  or  at 
least  was  much  smaller  than  the  "k  " -concentration  effect  itself. 
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In  line  with  this  finding,  Riedel  (39)  further  proposed  (in  the  study 
previously  cited)  that  the  same  assumption  should  hold  true  for  any 
binary  electrolytes  considered.  Thus,  the  reduced  thermal  conductivity 
correlation  for  aqueous  NaCl  solutions  shown  in  Figure  4.8  might  be  pre- 
sumed  to  hold  true  over  a  wide  range  of  temperatures,  rather  than  merely 
at  20°C . 

(3)  Pressure  Effect 

Mo  data  apparently  are  available  on  the  thermal  con¬ 
ductivities  of  aqueous  electrolytes  at  superatmospheric  pressures.  In 
view  of  the  relatively  weak  effects  of  concentration  and  temperature  on 
reduced  thermal  conductivity  "kr"  (as  compared  to  the  corresponding 
effects  on  reduced  viscosity  nur") ,  it  might  be  reasonable  to  assume  that 
the  presence  of  a  dissolved  salt  in  water  leads  to  no  special  effect  of 
pressure  on  thermal  conductivity,  and  that  pressure  dependence  for 
aqueous  electrolytes  can  be  satisfactorily  approximated  by  that  for 
pure  water. 

The  variation  in  the  absolute  thermal  conductivity  of  pure  water 
with  temperature  and  pressure,  as  reported  by  Bridgman  (4),  is  shown  in 
Figure  4.9.  It  is  seen  that  the  thermal  conductivity  of  water  increases 
only  by  about  20  percent  as  temperature  is  raised  from  0°C  to  100°C  (and, 
in  fact,  goes  through  a  maximum  in  the  temperature  region  above  100°C) . 

It  is  further  seen  that  thermal  conductivity,  at  various  temperatures, 
increases  by  about  6  percent  with  an  elevation  in  pressure  to  1000  atm. 

The  pressure  effect  is  roughly  comparable  to  that  observed  in  the  case 
of  viscosity,  but  the  influence  of  temperature,  obviously,  is  much  less. 

b.  Thermal  Conductivity  Correlations 

(1)  Relations  of  the  "Jones -Dole"  Type 

If  the  temperature  and  pressure  dependence  of  the  thermal 
conductivity  of  aqueous  electrolytes  are  taken  to  be  that  of  pure  water, 
the  only  remaining  correlation  of  interest  is  that  which  relates  thermal 
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FIGURE  4.9.  EFFECT  OF  TEMPERATURE  AND  PRESSURE  ON  THERMAL  CONDUCTIVITY  OF  PURE  WATER 


conductivity  to  concentration.  Presumably  a  semi -empirical  relation  of 
the  "Jones -Dole"  type  might  apply: 

kr  Z  1  +  a'Vc  +  b'C  +  d’C2  (4.25) 

in  which  kf  is  reduced  thermal  conductivity, 

C  is  concentration,  moles/liter, 

a',  b ' ,  and  c'  are  coefficients,  having  the  units 

of  (moles/liter) (moles/liter) ,  and  (moles/liter)-2, 

respectively. 

Inasmuch  as  thermal  conductivity  for  most  aqueous  electrolytes  decreases 
with  rising  concentration,  it  should  be  expected  that  one  or  more  of  the 
coefficients  a',  b' ,  and  d'  will  prove  to  be  negative  in  value. 

A  simplified  version  of  equation  4.25  considered  by  Riedel  (39)was 
based  upon  the  finding  that,  for  most  neutral  salts,  the  decrease  in 
thermal  conductivity  is  very  nearly  linear  in  molar  concentration: 

kf  -  1  +  b'C  (4.26) 

where  b'  is  understood  to  be  negative  for  most  electrolytes.  Riedel  then 

split  the  b'  coefficient  up  into  b'  ^  contributions,  based  upon  b'^a+  -  0, 

in  a  manner  analogous  to  that  previously  discussed  for  the  B  coefficient 

in  the  Jones-Dole  equation  for  viscosity.  This  assumption  of  the  additivity 

of  b'  contributions,  at  least  for  binary  electrolytes,  was  deemed  to  be 
ion 

valid  on  the  basis  of  the  relatively  extensive  amount  of  data  available 
for  such  systems. 

As  was  done  in  the  case  of  viscosity,  the  assumption  of  additivity 
of  individual  b'^on  contributions  was  applied  to  the  sea  water  system.  The 
resulting  prediction  is  shown  in  Figure  4.8 for  comparison  with  the  experi- 
mental  NaCl  values  and  the  previous  estimation  of  Krummel  (26) .  For  a 
concentration  of  3.5  weight  percent  salts  in  solution  (S  -  35  gm/kg) ,  the 
predicted  reduced  thermal  conductivity  of  sea  water  is  0.9937,  compared  to 
the  observed  value  of  0.9940  for  aqueous  sodium  chloride  at  the  same  concen¬ 
tration.  In  contrast  to  this  close  correspondence,  the  "kr"  value  estimated 
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for  sea  water  by  Krummel,  based  upon  the  assumption  of  equal  thermal 

dif fusivities ,  was~0.9575.  Apparently  this  crude  an  assumption  is  invalid 

for  this  type  of  comparison. 

Litvinenko  and  R' Jchenko  (29)  also  tested  an  equation  of  the 
"Jones -Dole"  type: 

kr  ^  1  +  a'Vc  +  b'C  (4.27) 


They  found  that  values  of  b  1  were  negative  for  nearly  all  aqueous  electro¬ 
lytes,  whereas  values  of  a1  were  small,  but  positive  in  a  few  cases.  Their 
results  will  be  further  discussed  in  the  section  to  follow. 

(2)  Correlations  Based  Upon  Apparent  Molal  Thermal  Conductivity 

A  second  method  of  relating  the  change  in  thermal  conductivity 
of  an  electrolyte  to  the  proportion  of  dissolved  salt  has  been  that  based 
upon  the  concept  of  apparent  molal  thermal  conductivity  defined  as  follows: 


®k 


(1000  t  iWk  •  lOOOkw 
m 


(4.28) 


in  which  0^  is  apparent  molal  thermal  conductivity,  cal/cm-aec-°C, 
m  is  molality,  moles  salt/1000  gms  water, 

M  is  molecular  weight  of  salt,  gm/mole,  and 
kw  is  thermal  conductivity  of  pure  water,  cal/cm-sec-°C. 


Kapustlnskiy  and  Ruzavin  (22)  attempted  to  show  a  linear  correlation  be¬ 
tween  0^  and  Vm,  in  which  values  of  the  extrapolated  limit, 0°, were  taken 
to  be  an  additive  function  of  the  contributions  of  the  ions  present. 
Prudnikov  (36),  on  the  other  hand,  demonstrated  that  this  type  of  corre¬ 
lation,  at  least  for  molalities  in  the  vicinity  of  1.0m  and  greater,  was 
inconsistent  with  the  experimentally-observed  linear  variation  of  total 
thermal  conductivity  with  concentration. 

Litvinenko  and  Radchenko  (29),  as  cited  previously,  proposed  a 
semi -empirical  relation  of  the  "Jones-Dole"  type  (equation  4.27)  between 
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electrolyte  thermal  conductivity  and  concentration.  They  also  examined 
the  "partial  molal  thermal  conductivity"  concept,  and  found  that  by  far 
the  greater  fraction  of  the  limiting  value  0^  could  be  ascribed  merely  to 
volume  changes  between  pure  water  and  solution,  and  with  thus  only  a  small 
portion  of  0°  was  actually  related  to  ion-solvent  interaction.  A  correspondence 
could  be  shown  between  this  latter  portion  of  and  their  empirical  b' 
coefficient  in  equation  4.27. 

In  view  of  the  difficulty  of  defining  a  total  mole  quantity 
in  a  multicomponent  mixture  such  as  sea  water,  the  apparent  molal  thermal 
conductivity  approach  appears  to  be  of  questionable  value  for  predicting 
the  behavior  of  such  systems. 

(3)  Correlation  Between  Thermal  Conductivity  and  Sonic  Velocity 
An  entirely  different  approach  to  the  correlation  of  thermal 
conductivity  of  aqueous  electrolytes  has  been  proposed  by  Eigen  (9)  ,  who 
based  his  method  on  a  demonstrated  correspondence  between  thermal  con¬ 
ductivity  and  sonic  velocity.  The  thermal  conductivity  of  pure  water  is 
considered  to  be  the  sum  of  two  terms: 

kw  -  k1  4-  kA  (4.29) 

in  which  k^  is  the  measured  thermal  conductivity  of  pure  water, 
cal/cm-sec-°C, 

k^  is  the  (theoretically)  calculated  thermal  conductivity 

of  pure  water  as  a  non-associated  liquid,  cal/cm-sec-°C,  and 

k.  —  the  difference  between  these  two  --  is  the  thermal  con- 
A 

ductivity  anomaly  of  water  arising  from  association, 
cal/cm-sec-°C . 

In  the  case  of  an  aqueous  electrolyte,  both  contributions  to  the  thermal 
conductivity  are  modified: 

k  *  ki'  +  kA'  (4.30) 

in  which  the  primes  denote  values  for  the  electrolyte  solution. 
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The  first  term  (kj  or  k^)  in  equation  4.29  or  equation  4.30  can  be 
shown  to  be  essentially  proportional  only  to  the  sonic  velocity: 


(4.31) 


in  which  U  is  sonic  velocity  of  electrolyte,  m/sec,  and 
is  sonic  velocity  of  pure  water,  m/sec. 

Inasmuch  as  the  sonic  velocity  in  most  electrolyte  solutions  is  greater 
than  that  in  pure  water,  k^'  may  be  expected  to  exceed  k^  for  pure  water. 

In  contrast,  the  second  term  (k^  or  k^' )  in  equation  4.29  or  equation  4.30  is 
decreased  in  direct  proportion  to  increasing  electrolyte  concentration: 


kA'  -  kA(l  -  Q'C)  (4.32) 

in  which  or'  is  a  constant,  (moles/liter) ,  and  C  is  concentration, 
moles/liter.  Values  of  the  constant  a'  can  be  calculated  for  simple  binary 
salt  solutions  if  specific  heat  data  are  available  for  the  electrolyte. 

As  an  approximate  application  of  this  correlation,  the  thermal 
conductivity  of  ~3.5  percent  (S  =  35  gm/kg)  sea  water  at  20°C  and  atmos¬ 
pheric  pressure  was  predicted  from  equations  4.29  through  4.32,  based  upon 
sonic  velocity  data  for  sea  water  as  summarized  by  Wilson  (47),  but  with 
the  constant  a  '  in  equation  4.32  calculated  for  binary  NaCl  electrolyte 
rather  than  for  sea  water.  The  result  of  this  prediction  was  an  estimated 
reduced  thermal  conductivity  (k^)  value  of  0.9928  at  these  conditions, 
which  may  be  compared  in  Figure  4.8  to  the  value  of  0,9937  predicted  from 
ionic  additivities.  Since  the  former  of  these  two  values  was  arrived  at 
by  taking  the  difference  between  the  positive  "k^"  correction  and  a  larger 
and  negative  "kA"  correction  (with  approximations  in  both  of  these),  the 
latter  "kr"  value  predicted  by  the  relatively  simpler  (albeit  more  empirical) 
method  is  considered  the  better  of  the  two. 


(4)  Sample  Predictions  of  Sea  Water  Thermal  Conductivity 

Inasmuch  as  completely  theoretical  methods  of  predicting 
thermal  conductivity  for  multicomponent  salt  solutions  appear  to  offer  no 


advantage  over  empirical  procedures  based  primarily  upon  experimental  data 
for  binary  salt  solutions  and  for  pure  water,  the  latter  are  recommended. 
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(a)  Predict  the  thermal  conductivity  of  ordinary  sea 
water  (S  35  gm/kg)  at  5°0  and  500  atm. 

k  of  H20  at  5°C,  1  atm  =  0.001428  cal/cm-sec-°C  [Figure  4.9] 

k 

of  H20  at  500  atm  *  1.031  [Figure  4.9] 

o 

k^  of  3.5  percent  sea  water  (20°C,  1  atm)  «=  0.994  [Figure  4.8] 

k  of  3.5  percent  sea  water  at  5°C  and  500  atm 

~  (0.001428) (1.031)  (0.994)  =  0.001463  cal/cm-sec-°C 

(b)  Predict  the  thermal  conductivity  of  10  percent  sea 
water  brine  at  60°C  and  1  atm. 

k  of  H20  at  60°C,  1  atm  =  0.001559  cal/cm-sec-°C 

kf  of  10  percent  sea  water  (20°C,  1  atm)  ~  0.982 

k  of  10  percent  sea  water  brine  at  60°C  and  1  atm 
«  (0.001559) (0.982)  -  0.001531  cal/cm-sec-°C 

4.5  THEORFTICAL  RELATIONSHIPS  BETWEEN  TRANSPORT  PROPERTIES 

Considerable  work  has  been  reported  in  the  investigation  of 
theoretical  relationships  between  viscosity,  electrical  conductivity,  and 
diffusion  coefficients  in  aqueous  solutions  of  electrolytes.  Thus  far,  those 
relations  which  have  been  developed  either  hold  true  only  in  extremely  dilute 
solution,  or  else  contain  empirical  elements  requiring  some  experimental 
measurements  before  predictions  can  be  extrapolated.  Thus  far  no  theoretical 
relations  have  been  suggested  between  electrolyte  thermal  conductivity  and 
any  of  these  other  three  transport  properties. 

The  simplest  type  of  theoretical  relationships  between  viscosity, 
electrical  conductivity,  and  diffusion  coefficient  are  those  based  originally 


[Figure  4.9] 

[Extrapolated 
from  Figure  4.8] 
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upon  Walden's  (46)  rule,  and  more  recently  outlined  by  Robinson  and 
Stokes  (40),  Harned  and  Owen  (16),  McLaughlin  (31)  and  others: 


(A. 33) 


D  a  A 


(4.34) 


in  which  ii  is  viscosity,  cp, 

A  is  equivalent  conductance, 


mhos -cm 


2  equivalent 
D  is  diffusion  coefficient,  cm  /sec,  and 


"or" indicates  proportionality,  but  by  an  unspecified  factor, 
and  where  A  approximately  can  be  defined  as: 


A  =  (4.35) 

in  which  L  is  the  specific  electrical  conductivity  of  the  solution,  mhos/cm., 
L°  is  specific  electrical  conductivity  of  pure  water,  mhos /cm. ,  and 
C*  is  electrolyte  concentration  in  equivs . /liter . 


Equations  4.33  and  4.34  can,  of  course  be  combined  to  yield  a  reciprocal  re¬ 
lation  between  viscosity  and  diffusion  coefficient: 


M.  or  (4.36) 

In  using  equations  4.33,  4.34,  or  4.36,  of  course,  it  must  be  borne  in  mind  that' 
these  relations  apply  only  to  very  dilute  (<0.01M)  solutions  of  electrolytes, 
and  for  a  given  temperature  and  pressure. 

For  electrolytes  of  higher  concentration,  both  equivalent  con¬ 
ductance  and  diffusion  coefficient  have  been  correlated  versus  concentration 
by  semi -empirical  relationships  of  the  "Jones-Dole"  type  (16): 

A  -  A°[l  +  AftVc  +  Bac]  (4.37) 

D  »  D°[l  +  ApVc  .  BjjC  ]  (4.38)* 

*  The  third  term  was  added  in  order  to  maintain  consistency  with  Eq.  4.37. 
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V 


in  which 


is  lim  A» 
C-*  0 


mhos -cm' 
equiv. 


D°  is  lim  D,  cm^/sec. 

C-*0 


A  ,  B. ,  A^,  B_  are  coefficients  having  appropriate  units. 

A  A  D  D 

In  line  with  the  limiting  rules  discussed  previously,  "A^"  and  "A^"  are 
both  negative,  and  both  have  essentially  the  same  value  as  "(-A)"  in  the 
conventional  Jones-Dole  expression  for  viscosity. 

The  similarity  between  these  semiempirical  laws  for  viscosity 
behavior  and  electrical  conductance  behavior  has  led  to  several  proposed 
relationships  between  the  two,  valid  over  wide  ranges  of  concentration 
and  temperature.  Falkenhagen  and  Leist  (11,12),  for  example,  have  shown 
that  calculable  viscosity  "corrections"  must  be  applied  when  predicting 
electrical  conductivity  for  electrolyte  solutions  in  the  region  of  high 
concentration.  Other  relationships  between  viscosity  and  electrical 
conductance  have  been  proposed  by  Kelbg  (23),  and  by  Izmailov  (17)  (the 
latter,  an  empirical  correlation  of  the  product  of  the  two  versus 
temperature) . 

A  somewhat  different  approach  has  been  taken  by  workers  in 
the  field  of  irreversible  thermodynamics,  including  Laity  (28),  and  more 
recently.  Miller  (32)  and  Lorenz  (30).  Theoretical  relationships  have 
been  derived  between  the  Onsager  interionic  interaction  coefficient 
L^  (“  1*21^  and  both  electrical  conductance  and  diffusion  coefficient. 

(The  L^2  coefficient  arises  when  transport  equations  are  written  in  terms 
of  total  ionic  "flows"  in  a  solvent-fixed  reference  frame,  taking  the  place 
of  the  more  conventional  individual  equations  for  mass  diffusion  and  elec¬ 
trical  conductance.)  Miller  (32)  has  pointed  out,  however,  that  no  such 
theoretical  relationships  can  exist  between  the  solution  viscosity  and 
either  of  these  other  two  transport  properties,  inasmuch  as  viscous 
phenomena  differ  from  diffusion  and  thermal  conduction  phenomena  in  tensorial  order. 
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The  conclusion  is  drawn  that  any  existing  empirical  correlations  between 
viscosity  (^)  and  either  D  or  A  for  electrolyte  solutions  of  appreciable 
concentration  must  be  strictly  accidental. 

Based  upon  all  preceding  considerations,  it  is  concluded  that 
no  suitable  completely  theoretical  relationships  presently  exist  between 
either  viscosity  or  thermal  conductivity  and  the  mass  diffusion  properties 
considered  in  Sections  2  and  3.  If  this  is  true  in  the  relatively  simple 
case  of  binary  electrolytes,  it  should  also  be  true  in  the  case  of  multi- 
component  mixtures  such  as  sea  water. 

4.6  SUMMARY  AND  CONCLUSIONS 

The  molecular  viscosity  and  thermal  conductivity  of  binary  and 
multicomponent  salt  solutions  have  been  investigated  by  recourse  to  the 
pertinent  literature.  Both  experimental  data  and  methods  of  theoretical 
or  semi-empirical  prediction  have  been  included  in  this  study.  The  primary 
application  of  these  data  and  methods  has  been  directed  toward  estimation 
or  calculation  of  the  transport  properties  of  sea  water  systems  at  various 
concentrations,  temperatures,  and  pressures.  A  secondary  emphasis  has 
been  given  to  methods  for  relating  these  two  bulk  transport  properties 
to  the  individual  molecular  ordinary  and  Soret  diffusion  coefficients. 

The  literature  consulted  for  this  study  included  more  than  100  references 
covering  work  mostly  reported  during  the  past  20  years. 

The  following  conclusions  may  be  drawn  from  the  studies  under¬ 
taken  to  date: 

(1)  The  viscosity  behavior  of  reasonably  dilute  aqueous  electro¬ 
lyte  solutions  greatly  resembles  that  of  pure  water,  in  that  temperature 
has  a  much  greater  effect  than  either  pressure  or  concentration.  Viscosity, 
in  other  words,  drops  by  a  factor  of  nearly  10  as  temperature  is  raised 
from  0°C  to  100°C,  whereas  pressure  effects  over  the  range  1-1000  atm 
generally  amount  to  less  than  ±  8  percent.  The  viscosity  of  a  3.5  weight 
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percent  solution  of  either  sodium  chloride  or  sea  water  Is  higher  than 
that  of  pure  water  by  an  amount  ranging  from  3  percent  to  8  percent, 
depending  on  the  temperature. 

(2)  The  thermal  conductivity  of  aqueous  electrolytes  generally 
appears  to  be  less  sensitive  to  either  temperature,  pressure,  or  concen¬ 
tration  than  viscosity.  In  the  case  of  pure  water,  thermal  con¬ 
ductivity  Increases  by  about  20  percent  as  temperature  Is  raised  from 
0°C  to  100°C,  and  Increases  by  about  6  percent  as  pressure  Is  elevated 
from  1  atm  up  to  1000  atm.  The  thermal  conductivity  of  aqueous  sodium 
chloride  and  most  other  similar  salt  solutions  Is  less  than  that  of  pure 
water,  but  the  relative  decrease  for  a  3.5  weight  percent  concentration 
is  less  than  1  percent . 

(3)  For  most  practical  purposes,  sea  water  can  be  considered 
to  be  essentially  a  3.5  weight  percent  solution  of  sodium  chloride,  but 
with  a  relatively  minor  proportion  of  magnesium  chloride.  In  ordinary 
engineering  calculations  the  bulk  transport  properties  (that  Is,  viscosity 
and  thermal  conductivity)  of  sea  water  can  be  estimated  reasonably  well 
from  available  data  on  sodium  chloride  solutions,  although  some  corrections 
based  upon  data  as  yet  unmeasured  might  be  necessary  In  the  case  of  more 
highly  concentrated  brines  where  the  solubility  limits  of  certain  minor 
salt  constituents  of  sea  water  would  be  exceeded.  (An  indication  of  available 
experimentally  measured  transport  properties  of  pure  water,  aqueous  sodium 
chloride,  and  sea  water  is  shown  in  Table  4.3)* 

(4)  Theoretical  relationships  do  exist  between  viscosity, 
electrical  conductance,  and  dlffusivlty  but  are  strictly  applicable  only 
for  the  case  of  very  dilute  electrolyte  solutions  in  which  the  ions  can 
be  considered  to  act  as  point  charges  which  are  effectively  not  hydrated 

by  the  solvent.  Any  relationships  between  viscosity  or  thermal  conductivity 
and  ordinary  or  Soret  mass  diffusion  coefficients  for  salt  solutions  In  a 
concentration  region  of  practical  interest  would  no  doubt  have  to  contain 
some  elements  of  empiricism. 
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TABLE  4.3.  AVAILABILITY  OF  TRANSPORT  PROPERTY  DATA 
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It  should  be  noted  that  the  results  reported  above  include 
the  first  attempt  to  apply  the  semiempir ical  Galinker  method  to  the 
prediction  of  the  viscosity  of  sea  water.  In  addition,  several 
techniques  were  used  to  predict  the  thermal  conductivity  of  sea  water 
and  the  usual  prediction  method  (equal  thermal  dif f usivities)  was 
shown  to  be  in  considerable  error. 
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SECTION  5 


GENERAL  SUMMARY  AND  CONCLUSIONS 


% 


1 


This  investigation  has  resulted  in  the  development  of  special 
high  pressure  cells  capable  of  measuring  the  ordinary  and  thermal  diffusion 
properties  of  binary  salt  solutions  for  pressures  in  the  range  of  1  to 
1000  bars.  The  equipment  was  used  to  measure  these  properties  for  the 
3.5  percent  by  weight  sodium  chloride— water  system.  The  results  were 
correlated  with  available  approximate  prediction  techniques.  All  binary 
ordinary  and  thermal  diffusion  coefficient  data  available  for  the  most 
predominate  salts  of  the  sea  water  system  were  compiled  and  analyzed. 

Studies  were  made  of  the  apparatus  required  in  order  to  make 
meaningful  measurements  of  the  ordinary  and  thermal  diffusion  properties 
of  multicomponent  aqueous  salt  solutions  such  as  sea  water.  Appropriate 
designs  and  operating  techniques  were  developed  based  upon  the  results  of 
these  studies. 

An  extensive  review  and  analysis  was  made  of  all  available 
literature  pertaining  to  the  viscosity  and  thermal  conductivity  of  salt 
solutions,  including  seawater.  Several  theoretical  prediction  methods 
were  applied  to  the  sea  water  system. 

On  the  basis  of  the  results  obtained,  the  following  general 
conclusions  may  be  made  (see  Sections  2,  3,  and  4  for  more  detail  and 
supporting  data): 
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(i)  the  effect  of  pressure  on  both  the  ordirtfcry«*dif fusion  coefficient 
and  the  Soret  coefficient  of  the  NaCl-l^O  system  at  25°C  in  the 
pressure  range  of  1  to  1000  bars  is  sufficiently  small  that 
it  can  be  neglected;  the  ordinary  diffusion  results  are  con¬ 
sistent  with  estimates  based  on  various  extrapolation  methods; 

(ii)  the  binary  diffusion  properties  of  the  major  salt  constituents 
of  sea  water  are  sufficiently  different  that  the  diffusion 
behavior  of  sea  water  probably  can  not  be  predicted  on  the  basis 
of  only  the  NaCl-l^O  system  behavior;  thus,  a  significant 
pressure  effect  could  still  exist  for  the  sea  water  system  even 
though  none  was  found  for  the  NaCl-l^O  system  (at  25°C); 

(iii)  the  thermal  diffusion  effect  in  the  NaCl-H^O  system  is  so 

temperature  dependant  that  at  about  7°C  the  Soret  coefficient 
would  be  expected  to  be  zero  (no  thermal  diffusion  effect)  and 
at  lower  temperatures,  the  Soret  coefficient  takes  on  negative 
values  —  that  is,  the  thermal  diffusion  effect  is  reversed  and 
molecular  salt  transport  can  now  occur  from  warm  regions  in  the 
fluid  to  the  colder  regions  instead  of  vice  versa. 

(iv)  the  viscosity  and  thermal  conductivity  of  dilute  aqueous  electro¬ 
lyte  solutions,  including  sea  water,  are  very  similar  to  the 
corresponding  properties  of  pure  water;  the  thermal  conductivity 
appears  to  be  less  sensitive  to  temperature,  pressure,  and  con¬ 
centration  than  the  viscosity;  engineering  estimates  of  these 
two  bulk  transport  properties  can  be  made  by  assuming  that  sea 
water  is  essentially  a  three  and  one  half  weight  percent  sodium 
chloride  in  water  solution. 

TVfo  major  areas  will  require  further  investigation.  These  are 
the  multicomponent  diffusion  problem  (oriented  towards  the  sea  water  case) 
and  the  problem  of  how  temperature  affects  the  diffusion  properties  at 
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temperatures  where  the  density  of  water  exhibits  a  maximum  with  temp¬ 
erature  (~4°C).  Further  binary  experiments  must  be  made  with  the 
NaCl-i^O  system  as  well  as  the  other  binary  aqueous  salt  solutions  which 
constitute  sea  water  in  order  to  confirm  the  results  so  far  obtained. 
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APPENDIX  A. 1 


PROCEDURE  FOR  ORDINARY  DIFFUSION  COEFFICIENT  MEASUREMENT 

AT  HIGH  PRESSURE 


1.  Disassemble  the  pressure  cell  completely.  Wash  all  parts  in  hot 
water.  (Use  wire  brush  to  remove  rust  before  washing).  Rinse  in 
distilled  water.  Dry  parts  by  acetone  rinse  and  air  dry. 

2.  Wash  diffusion  diaphragm  in  Special  Vacuum  Washing  Apparatus,  using 
distilled  water.  Rinse  in  acetone,  and  dry  by  pulling  air  through. 

3.  Boil  stock  NaCl  solution  under  vacuum,  and  a  supply  of  distilled  water 
under  vacuum. 

4.  Check  condition  of  diaphragm.  Measure  bomb  temperature.  Place  dia¬ 
phragm  in  beaker  of  degassed  NaCl  solution.  Put  beaker  under  bell 
jar,  and  evacuate.  Allow  to  stand  under  vacuum  for  1  hour  -  shaking 
occasionally  to  dislodge  bubbles.  Flush  diaphragm. 

5.  Assemble  pressure  cell  in  following  order: 

(a)  install  bottom  closure  piece,  packing  washers,  etc.  first. 
Place  temporary  plug  in  pressure  inlet  -  making  sure  the 
capillary  section  is  in  place  in  the  pressure  inlet  hole. 

(b)  Place  cell  on  stand  and  fill  bottom  section  with  NaCl 
solution  above  diaphragm  seat. 

(c)  Loosen  temporary  plug  on  bottom,  and  allow  a  few  drops  of 
salt  solution  to  run  out.  Retighten,  and  note  gas  bubble 
evolution  in  cell.  Repeat  until  no  more  gas  is  evolved. 
Dislodge  any  air  bubbles  from  the  cell  side  wall. 
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(d)  Press  diaphragm  into  place,  through  the  NaCl  solution, 
using  a  piece  of  plastic  pipe. 

(e)  Using  a  plastic  syringe,  remove  excess  '’aCl  solution  from 
above  the  diaphragm.  Be  careful  not  to  touch  diaphragm 
with  the  syringe.  Rinse  top  cell  section  4  times  with 
distilled  water,  removing  each  rinse  before  adding  the 
next.  Do  not  pour  rinse  directly  on  the  diaphragm. 

(f)  Fill  the  tcp  half  of  cell  with  distilled  1^0  to  inch 
above  the  sealing  step.  Dislodge  bubbles,  insert  top 
sleeve  liner,  and  squirt  water  between  it  and  cell  wall 
to  dislodge  bubbles. 

(g)  Start  the  run  timer. 

(h)  Assemble  the  top  closure  piece,  packing  washers,  drive  nut, 
and  knurled  retaining  cap.  Insert  smoothly,  and  screw  down 
the  driving  nut.  (Remove  excess  water  as  it  is  forced  out 
the  pressurizing  line.) 

(i)  Insert  the  capillary  section  in  top  pressurizing  hole. 

6.  Move  assembled  cell  to  bath,  and  connect  pressurizing  lines  to  top  and 
bottom  pressure  inlets.  Tape  lines  to  cell  body  for  alig  ^ent. 

7.  Place  cell  in  bath  and  start  stirrer. 

8.  Pressurize  cell  over  desired  time  interval  with  the  displacement  pump. 

No  stirring  until  pressure  is  reached. 

9.  Allow  stirring  to  continue  for  the  desired  time  measuring  the  electrode 
pair  resistances  at  regular  time  intervals.  Depressurize  over  the 
desired  time  period. 

10.  Stop  stirrer  and  remove  ceil  from  the  bath. 

11.  Disconnect  pressure  connection  at  top  of  cell,  and  insert  temporary 
plug.  Turn  cell  over,  and  remove  bottom  pressure  connection. 

12.  Install  retainer  ring  and  remove  bottom  closure  piece,  washers,  etc. 

13.  Using  a  plastic  syringe,  remove  bottom  solution  and  transfer  to  sample 
bottle. 

14.  Invert  the  cell  and  repeat  with  top  section. 

15.  Disassemble  cell,  and  wash  all  parts  thoroughly.  Rinse  with  distilled  water 

and  acetone.  Air  dry.  Check  condition  of  diaphragm. 
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APPENDIX  A. 2 


APPROXIMATE  ESTIMATE  OF  THE  BULK  TRANSPORT  OF  SOLUTION 
IN  THE  ORDINARY  DIFFUSION  BOMB  ON  DEPRESSURING 


For  Experiment  ODBRS-13*  (t^  to  t^,): 


Diaphragm 


Volume  fraction  shoved 
into  top  chamber  =  a 


Mole  Balance  for  Top  Chamber: 


aV 


B 


V  +  [VT  -  “  VB  ]  CT,  '  VT  CT„ 


(A.  1) 


Mole  Balance  for  Bottom  Chamber: 


(1  -  or)  VB  C,  -  VB  CBfi 


But  VT  m  Vg,  therefore, 


cr  C«  +  (1  -  cr)  Qy  -  C, 


B. 


f  Tf 


(1  -  at)  C.  -  C. 

Bf  Bf, 


(A.  2) 


(A. 3) 
(A. 4) 
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The 


These  two  equations  may  be  solved  for  the  unknowns, 
quantity,  C  ,  is  then  found  from  (see  Eq.  14), 

D 

O 


o  and  C 

B 


f 


(A. 5) 


For  Experiment  ODBRS-13**  (t^  to  t^,): 


Volume  fraction  shoved 
into  bottom  chamber  =  o 

Diaphragm 


Mole  Balance  for  Bottom  Chamber: 


a  V*"  €<"  +^B-oVT]C  -  VRCB 


T  T, 


B, 


f  1 


(A. 6) 


Mole  Balance  for  Top  Chamber: 


(1  -  a)  VT  CT  -  V„  C 


T  T 


f 


But  VT  cs  Vg,  therefore, 


a  CT  +  (1  -  or)  C.  »  C 


Bf  Bf  i 


(1  -  a)  C_  *  C_, 

Tf  V 


(A.  7) 


(A. 8) 


(A. 9) 
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These  last  two  equations  may  be  solved  for  the  unknowns,  a  and  C_  . 

-  Bf 

C_  is  found  from  Eq.  A. 5.  D19  can  then  be  calculated,  using  the 


electrode  determined  values  of  C_  and  C_  ,  from  Eq.  12. 

It  should  be  noted  that  C_  an  C_  are  the  usual  pipette 

Tf '  Bf ' 

cell  measured  concentrations  measured  for  the  solutions  withdrawn 


from  the  ordinary  diffusion  bomb  after  depressurizing  to  1  atm. 
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NOMENCLATURE 


°Fick 


Er  (  > 


coefficient  in  Eq.  4.11  (moles/liter)  ^ 
height  of  binary  Sortt  cell  (cm) 
coefficient  in  Eq.  4.25  (moles/ liter)  ^ 

diaphragm  pore  cross-sectional  area  available  for  diffusion 
(cm  ) 

"Jones-Dole"  coefficient  (moles/liter)  * 

coefficient  in  Eq.  4.11  (moles/liter)  * 

coefficient  in  Eq.  4.25  (moles/liter)  * 

"Jones-Dole "coefficient  (moles/liter)  * 

property  parameter  in  Eq.  18 

coefficient  in  Eq.  4.9  (moles/ liter)  * 

coefficient  in  Eq.  4.11  (moles/liter) 

molar  concentration  (moles/liter) 

equivalent  concentration  (equivs/1 iter) 

„2 

coefficient  in  Eq.  4.25  (moles/liter) 

_2 

"Jones-Dole"  coefficient  (moles/liter) 

.  2 

coefficient  in  Eq.  4.9  (moles/ 1 iter) 

2 

Fick  diffusion  coefficient  (cm  /sec) 

2 

binary  ordinary  diffusion  coefficient  (D^  *  (cm  /®ec) 

function  defined  in  Eq.  4.15  (Kcal/mole) 
activation  energy  for  viscous  flow  (Kcal/mole) 

Faraday's  constant 

a  function  of  C_  and  C_  (see  Eq .  16) 

T  c 
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2  2 

■  flux  density  of  species  2  (tnole/ctn  sec,  or  gtn/ctn  sec) 

■  thermal  conductivity  (cal/cm  sec°C) 

-  thermal  conductivity  anomaly  of  pure  water  [Eqs.  4.29 
and  4.32]  (cal/cm  sec°C) 

•  thermal  conductivity  anomaly  of  electrolyte  [Eqs.  4.30 
and  4.32]  (cal/cm  sec°C) 

■  ideal  thermal  conductivity  of  pure  water  [Eqs.  4.29 
and  4.31]  (cal/cm  sec°C) 

•  ideal  thermal  conductivity  of  electrolyte  [Eqs.  4.30 
and  4.31]  (cal/cm  sec  °C) 

■  constant  in  Eq.  4.14  (cp) 

•  function  defined  in  Eq.  4.15  (cp) 

•  effective  diffusion  path  length  through  diaphragm  (cm) 

•  specific  electrical  conductivity  (mhos/cm) 

-  Onsager  interaction  coefficient  (units  depend  on  associated 
driving  force  and  flux) 

■  molality  (moles/1000  gms  1^0) 

■  molecular  weight  (gm/gm  mole) 

■  molarity  (moles/liter) 

•  mole  flux  of  species  2  relative  to  stat ionary  coordinates 

2 

(mole/cm  sec) 

23 

«  Avogadro's  number  (6.023  x  10  molecules/gm  mole) 

■  pressure  (atm  or  bar) 

■  gas  constant  ■  0.001987  (Kcal/mole  °K) 
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R  *  electrical  resistance  between  a  pair  of  electrodes  (ohm) 

*  radius  of  diffusing  solute  particle 

S  *  salinity  (gm/Kg  1^0) 

*  degree  of  saturation 

t  *  time  (seconds) 

T  ■  temperature  (°K) 

v^  ■  velocity  of  species  i  (cm/sec) 

U  *  sonic  velocity  (m/sec) 

x  ■  mole  fraction 

x  ■  distance  measured  up  from  cold  wall  in  binary  Soret  cell  (cm) 

X  *  constant  defined  in  Eqs.  4.6  and  4.7 

■  thickness  of  section  i  in  Soret  bomb 

Y  -  constant  defined  in  Eqs.  4.6  and  4.8  (moles/ liter)  * 

Y  *  weight  percent 

y  *  mass  fraction 

Zi  ■  ionic  valence  (including  sign)  of  ion  i  (z^  ■  0  for  electrically 

neutral  species,  such  as  l^O)  (  |z |  ■  number  of  protonic  charges 

per  ion  *  number  of  gm  equivalents  of  ion  i/gm  mole  of  ion  i) 


Greek  Symbols : 


constant  defined  in  Eq.  4.13 


or  *  net  volume  fraction  of  solution  in  ordinary  diffusion  chamber 

which  is  bulk  transported  across  cell  diaphragm  due  to  depressuring 
process  (see  Appendix  A. 2) 


or  •  signifies  proportionality  in  Eqs.  4.33  and  4.34 
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a' 

a 

a 


x 

e 


i 


A 


P 

TT 

* 

Y 


V 

CT 

a 


o 

i 

12 


constant  defined  in  F.q.  4.32  (moles/liter) 

-2 

cell  constant  (cm  ) 
constant  defined  in  Eq.  4.13 

difference  between  the  value  of  the  argument  at  pressure  P 
and  that  at  1  atm 

2 

equivalent  conductivity  of  ion  i  (mho  cm  /gm-equivalent  of  ion  i) 

characteristic  ordinary  diffusion  relaxation  time  (sec)  * 

2 


tt2D 


12 


equivalent  conductance  p 


mhos  cm 
gm  equivalent. 


0 


viscosity  (cp  or  gm/cm  sec) 

3 

density  (gm/cm  ) 

3.1416 


partial  molal  quantity  (units  vary  with  application) 
function  defined  in  Eqs.  4.10  and  4.11 
coefficient  defined  in  Eq.  4.11 
multicomponent  ionic  Soret  coefficient  (°K  *) 
binary  Soret  coefficient  (a^  "  "a21^  (°K  ^ 


Subscripts: 

B  ■  bottom  chamber 

D  •  denotes  coefficient  in  Eq.  4.38 

DP  *  diaphragm  pores 

F  ■  feed  solution 

f  *  end  of  ordinary  diffusion  experiment  before  pressure  has  been 

reduced 
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f' 


i 

ion 

k 

m 

o 


P 

r 

r 

s 

T 

w 

w 

A 

0 

1,2 


1 ,2  | . . . 


end  of  ordinary  diffusion  experiment  after  pressure  has  been 
reduced  to  1  atm 

species  i 

denotes  coefficient  for  individual  ion  contribution 
denotes  function  with  respect  to  thermal  conductivity 
Monel 

start  of  ordinary  or  thermal  diffusion  experiment  when  pressure 
in  bomb  has  been  established 

start  of  ordinary  diffusion  experiment  before  pressure  has  been 
increased  from  1  atm 

denotes  property  at  indicated  pressure  P 

denotes  "reduced"  property  (ratio,  electrolyte:  pure  water) 
phenolic  resin 

denotes  property  at  saturation 
top  chamber 

denotes  property  for  pure  water 
NaCl-l^O  solution 

denotes  coefficient  in  Eq.  4.37  for  equivalent  conductance 
denotes  property  at  zero  pressure 

denotes  two  different  concentrations  in  Eqs.  4.11  and  4.12  — 
more  generally,  these  subscripts  denote  species  1  and  2; 
species  2  is  taken  to  be  the  solute  species 

denotes  number  signifying  each  component  in  Eqs.  4.16  through  4.23 
stationary  state  (t-*®) 
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Superscripts; 


-  quantity  associated  with  the  beginning  and  end  of  an  ordinary 
diffusion  experiment  when  the  pressure  is  at  1  atm 

0  -  denotes  limiting  value  for  infinite  dilution 

m  ■  mass-average  reference  system 

N  ■  mole-average  reference  system 


Overbars : 

■*  ■  vector 

-  ■  mean  or  average;  integral  property  (Integrated  with  respect 

to  time  and  concentration) 


Mathematical: 

7  •  "del."  or  "nab la"  operator 

In  ■  natural  or  naperian  logarithm 
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